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1. In support cf the US^USSa exchange, we have 
prepared two documents. 

a. A geographical description of the 
southern Sierra Nevada/ San Joaguin Valley 
test site was sent to the NASA Technical 
Monitor in June 1977. 

b, A preliminary Gescription of the 
Soviet test site, the Fergana Valley, was 
sent tj the NASA Technical Monitor in April 
1977. y. In August 1977 we received a 
transliition of the aassian description of the 
Fergana area, and we have apdated our 
original descripticn tc include this new r? 
material. This versicn is included as 
chapter VI of this report. 

2. »6 have alsc made seme steps toward the 
development of an ehergy balance snowmelt model for 
large basins. This werk will be supported in the 
future by NASA Grant NSG-5062, NCAA, and the 
Oniversity of California Water Resources Center, 
Therefore we have deveted mest of our resources this 
winter to maintenance of the data coilectxcn 
programs, rather than tc development of the models, 

a, Ke have extended a surface climate 
simulation model over rugged terrain, and 
have tested it with soil temperature data 
from an area in the Wratgell Mountains, 

Alaska. These data were the most extensive 
we could find fer a mountainous area. A 
paper cn "An Approach toward Energy Balance 
Simulation over Rugged Terrain" has been^^^^ 
submitted for publication, and is included as 
Appendix 1 of this report. 

• i ' ' 

b . We have made progress toward the 
develcpment of a tepegraphic spectral solar 
cadiat icn model, especially designed for 
snowmelt iDcdelinyl These erforts are 
described in chapter II, and the data 
collected are in Appendix 2. 

o, Se have made progress toward 
determination of snow reflectance from 
Land sat data. These efforts are described in 






Ci 


d. Cue precress teward a topographic 
longwave radiatico model is described in 
chapter IV, and the data ccllected are in 
Appendix 3. i 


e. Cur investigations into snow-soil 
heat and aoisture flow in the southern Sierra 
are described in chapter V» Data collected 
toward this effort are in Appendix 4« 
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II. Developiisnt cf a lopoyraphic Solar 

% Hadiatioa Alyocltha roc Snowaelt Hodeling 

1 

' 

In high^ Douiitainous areas with sparse vegetation^ a 
najor source of energy for snow aetaeorphrsa ind aelt is 
inccsing solar radiation. Xo such environments melt can take 
place af lew temperatures, even kelcw freezing, and 
conventicnal index snowmelt iiiodels may perform iadly. 

aecause of the characteristics of snow and of 
mountainous terrain, a solar radiation algorithm should have 
the fcllcwinq attributesi 

1. It should be spectral, because of the 
^ variation of snow reflectance with wavelength. 

2.. It should account for radiation variatiens 
due to slope, exposure, and horizons. 

3. It should include calculations for radiation 
reflected from adjacent terrain. 

4. It should be able to make estimates of the 
reguired atmospbecic atteo'uation parameters from a 
limited set of measurements. 

5. It should previde for attenuation by clouds. 

6. It should be computationally economical. 

Thus far our werk has been entirely on a clear-sky 
version of a solar radiaticr atodel. For clear-s,k.y conaitiens 
we have sufficient data to test cur model, but we dc not have 
sufficient clcudy-sky data. The testing of our algorithm is 
still incomplete; therefore here cnly the basic structure is 
presented . , 

fte use the NCIC Digital lerrain Tapes to calculate 
slope, exposure, and horizon information for any point in a 
^|d4.gitiz6d grid. The method is described in detail in a paper 
by^^js^zier and Outcalt, which has been submitted for 
pl^bliciticn, and which is iacluded as an appendix to this 

(f^^port..X, 



A spectral solar radiation acdel modified, froiE Giergis 
(1975) is used to estimate beam and diffuse radiation. The 
variable atmospheric parameters are ozone, precipitable water 
vapor, and the Beta coefficient in the Angstrom aerosol 
scatterxnq function. Values for waveleny th-depende nt ozone 
absorption are taken from Inn and Tanaka ( 1953) and Leighton 
(1961). Water vapor absorption values are taken from Gates 
and flarrop ( 1963) . Non-vatiahle absorpticn coef f icidnts fot 
the red molecular oxygen bands (leighton, 19b1) and the 
near-infrared methane and carbon dioxide bands (Gates and 
Har rep, 1963) are also included. Va>l ues for Jjayleigh 
scattering are taken from Penndotf (1957). The exponent in 
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the Angstccm aecosol scattering functiou is assumed to be 
-1.3, although this assuopticn can be laodified if carrcw band 
spectral measureBents are available. 

Solar gecmetry is calculated by standard methods (e. g. 
Sellers, 1965), and path lengths tor air mass and water vapor 
are calculated by Kasten’s (1956) methods. The method for 
ozone path length calculation is from Lacis and Hansen 
(197b). We use the values tabulated by Makarova and 
Kharitinov (1S72) for the solar ccnstant. 

Diffuse and bacXscattered radiation are calculated by 
Giorqis* (1975) aethcdi A pcovisicn for specular reflection 
from adjacent slopes froB Paltridge and Platt (1976) is 
included. However, we have only approxamations for the 
specular ccmpcneht cf snow reflectance. Thus this is a 
source of possible error. 

We are new Bcdifying the pregram so that two 
measurements, taXen at : different solar zenith angles, can be 
used tio estimate the necessary ataospheric parameters by 
solving a set of non-linear algebraic eguatious using Ecown's 
(1969) method. An earlier approach, whereby we tried to 
define and measure integrated absorption and scattering 
coefficients over each Landsat band, proved infeasible, 
perhaps because the band widths are too large. Hence we have 
adopted the approach described above, which is independent of 
the particular wavelength interval of the measuring 
instrument, A remaining unsclved problem is how to ccrrect 
for systematic and random measurement errors. 

We will use out extensive solar radiation ffleasurements, 
corresponding to the four landsat visible and near— infrared 
bands, taXen during the 1977 and 1976 snow seasons, to test 
our model. These data are included in Appendix 2. 
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III. Bstiaating the Shcctwatie Albedo of aa Aljfirie 
Snow Covec Using Satellite Badiauce Measureaects 


Introducticn 


Knowledge of the spatial and temporal variation of 
shortwave albedo is essential for energy balance computations 
over an alpine snow cover. Ffeguent/ 'comprehensive field 
neasurements of surface reflectance over an area the size of 
the southern Sierra are impractical owing to the area's 
spatial diversity, extent and relative inaccessibility during 
the snoe season. Derivation of an acceptabie estimate of 
surface reflectance from satellite radxance aeasurements 
appears to be the nest cost'>e£t6ctlve neans of ccapleting the 
shortwave cadiaten inputs tc an euergy-balance snowmelt 
model. 


Relevant Characteristics cf Landsat liSS Data 

1. Spatial characteristics 

a. sanpliEg interval: The nominal IF07 of the 
Landsat MSS is 79 x 79 metars on the ground. along 
scan oversampling yields a mean sample spacing of 56 
meters in the x-dixecticn. Thus every 15-minute USGS 
guadrangle {the fundamental unit in cur geograhic data 
base) contains well in excess of 1Q5 Liadsat pixels # 
which poses a significant cempu tationai problem if each 
pixel is to be evaluated individually. On the ether 
hand, the sampling iDtexval may be sufficiently great 
that it Basics seme of the effects of terrain and 
veqetaticn cn reflectance. 

i : , : 

b. areal coverage: Three frames of Landsat 

imagery are reguired fer complete coverage ;bfi thei sfudy^- 
area. However, less than cDe-third qf any given frame 
covers an area cf interest. Extraction and 
concatenation of the study. area coverage from each 
frame into a single study area dataset is reguired, for 
each frame triplet. 

c. distorticn; landsat fiSS imagery inherently 
contains numerous geometric distertions iVan Sie anii 
Stein, 1976; Bsmstein and Fecneyhough, 1975) which 
ccmplicate the assignment of pixeirddiance values to 
specific Iccaticns on the gtcund. Some of these 
distortiens are systematic and may he cemovea by 
mcdelingf ethers regrirs inf ermation regarding 
spacecraft attitude, altitude, and position which is 
net directly available to a sufficient level of 
accuracy, and must therefore be statistically derived 
with respect tc g rcund central points in the image 
whose gecdetic locat-iciis are known. Ins desired result 
Of these two parallel, proceiures is a set of raapping 
trans.torms rela-i-lng amage and geograhic coerdinate 
systems./':. 
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2. : Radioaetric cbacactecistics 

a. dynamic range: laodsat dSS data frequently 

exhibits saturaticn over sacfc (ilcaianis et al, 1975). 
Such saturated pixels are esseatially laeless in 
radiance calculations. £y cut oi»n estimates 
(inspection cif thresbcld shade prints) , up to one-third 
o£ the most heavily sncwccvered areas may saturate in 
MSS bands 4 through 6* Little saturation seems to 
occur in Band 7, which is in accordance with labciatory 
measurements of the spectial reflectance of snow 
(C'Brien and Munis, 1S75). 

b. atmospheric effects: Any spaceborue sensor 

measures a combinaticn cf surface and atmospheric 
radiance.. Surface upuelllng radiance is subject to 
selective atmospheric attenuation. The atmosphere 
itself reflects back into space a portion of the 
incident sclar radiation. The magnitude of this 
combined attenuation and addition must be calculated 
for each Lacdsat MSS '<iavelength band. In addi-tion, 
some correction must be made for the change in path 
length as a functicn cf earth curvature and scan mirror 
angle. 

3. Spectral characteristics 

a. identification ox snew: Snow has a higher 

albedo than almost any cth|)r phenomenon in the Landsat 
ass bands, and thus may he readily identified, except 
in cases where snow and cloud cover lare coincident. It 
has been demonstrate d {Eattclucci et al, 1975) that 
snow is virtually indistingaishabie from clouds in the 
Landsat BSS bands. A means must be developed to 
accurately determine which pixels represent snew cover 
before radiance values are used in furxher 
computations, 

b, applicability tc snow albedc mapping; The 
Landsat MSS bands span that pcrtion of the spectfum = 
where the albedo of snow exhibits the greatest temporal 
variation (US Army Corps of Engineers, 1956) , and thus 
requires the most freguent updating for input to an 
energy-balance model. 

4. Temporal characteristics 

landsat's 9-day repea-t cycle (combinaticn cf 
landsats 2 and 2) is far frcit optimum for snowcever 
reflectance monitoririj. Since the most pronounced 
changes in albedo occur most nearly following a new 
snowfall, such a snowfall occurring between overpasses 
would be quite inaccurately represented, noth in terms 
of areal extent and albedo, by extrapolation from the 
previous overpass. 


Data Collectlcxi and fcccessioj 
1. RaH data 

a* satellite data: Since oost o£. ouc analysis 

techniques ace digital, we ch tain our satellite data as 
CCTs, tor those ovec^ass dates ^'pc which ue have 
c, ^incident ground truth aeasuceaiints. 


b. rield measucenientsi Cur shortwave ground 

truth consists oc a sequence of oeasuceaents at each 
site using an Sxctech 100-A landsat Gcoand Truth 
Radioneter. We have concentrated on the iollcwing four 
field test sites: Gotten \^ood Eass-Horseshoe aeadow 

(Kern river basin). Bullfrog Lake-7idette aeadew (Kings 
river basin) , Bishof tass (San Joaquin river basin) , 
and aaomoth Mountain (adjacent to San Joaquin rivet 
basin). The first three areas ate sufficiently resote 
tc preserve relatively undisturbed snow surfaces, yet 
ate accessible in che to two days on skis. The Mamnioth 
area is the location of ouch present and developing 
instrumentation by vaticus agencies in cooperation with 
tbe ski area managemeht. All four areas are easily 
located on landisat iiaacery and have measurement sites 
which are relatively icee of tcpograhic effects. 

c, topogriphic data: 5ie have developed software 

to extract 15-minate guaurangles from NCIG Digital 
Terrain Tapes, and resample them to any desired spatial 
rescluticn. The guadrangles from the funaanental grid 
of our study area database, and all Landsat and fiefd 
measurements are Iccated with respect to the terrain 
grid. * 


d. ground ccntrel data: Ground control points 

ate located on both USGS 15-minute guadrangle macs and 
on line printer shade prints of the Landsat imagery. 

«6 ace in the process of iitplementing a digital video 
display method of identifying control points in the 
image. 


2, Simulated data 


a. atmospheric parameters; Atmosphetic path 
radiance and transmissivity are calibrated with 
radiometer field measurefflents, then simulated by a 
model described alsewhere in this paper. 

b. solar irradiance at the sarface; Beam and 
diffuse solar radiaticn axe sittulted using a 
ccmhinaticn cf the at Bcspheric model, terrain data, and 
an ephemeris pregram which we have ueveioped. 


3. Data processing 

The principal data processing prooieni is the precise 
location of field ttasm-enient sites on the Lapdsat CCTs. This 
has been partially alleviated selecting test sites that 
are uniform and extensive enough so that a moderate 
locational error (eg, < 2 pixels) iill still jfield acceptable 
results. For areas where tcpcgrahic effects are important, 
however, the positic ral arro r must be less than one pixel. 

Se are currently perfcrming image rectification using 
the DIBS program suite deveicpea at GSFC. This has proven 
unsatisfactory for three reasons. First, the program runs in 
a batch mode, which leads to unacceptably long turnarcund 
times on a heavily-used educational system, especially when 
generating shade prints for GCP location, which reguires a 
special configuration of out system's Irne printer. Second, 
the programs are poorly accuaented, difficult to modify, and 
procuce cryptic output and diagnostics. Third, the programs 
ate collectively quite expensive on our system conf iguraticn 
(for example, reformatting and deriving napping functions 
costs about $200 in computer time alone), rfe have thus begun 
implementing a cut-down, highly specialized rectification 
pcogtam on a minicomputer (?DP-1 1/45) which has significantly 
more online disk storage (60C Mbytes) and an interactive 
operating system (OSIX) » 

A secondary data processinvj problem, which is celatively 
str aightf crwa rd , is the ccnversicn of the Landsat radiance 
numbers to sarfaceupwelling radiance, using a combination of 
a table lcckup„and the atBcspheric model. 


Initial Verification 

Although both Landsat and field data are available for 
several dates during 1976-77 snow season, a pteliminary 
calculation of albedo based on space measurements has not 
been attempted for the following reasons: 


The initial version of the atraospheric model 
d numerical instabilities whicn necessitated 
g tedesigneg, ^hich is still under way. 

2, Several erccrs were discovered in the field 
radiometer measurements, necessitating a statistical 
filtering of the recorded values. 

3. The image rectification prpgrams were 
initially run on a summer image, to facilitate 
Iccation of control points in the study area. We have 
had difficulty registering a winter image to the 
rectified summer imace. The programs we are using 
also do not readily rotate imagery so that it is 
aligned to the cardinal directions, as are our terrain 
datasets and USGS maps. 


1 . 
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fie feel that these eccocs have been ceuedied by aoving 
to a Dev computer envircQaenti developing a uev atncspheric 
aodel, and recalibrating cur field instruaents, 

dnce vs are satisified that we can compute snow albedo 
over a uniforn snowpack, we will attsapt to extend the model 
over a topcgrahlc surface, which nay necessitate corrections 
for the high specular reflectancs wlich snow exhibits at low 
sun ingles (Pirahlrn and Eaton, 1974), Areas where this is a 
problem may |fa€ readily identified by combining the epheneris 
calculations (prcgram) with slope and aspect inforsation 
derived from the terrain tapes. Existing numerical 
corrections for anisctrcpic reflectance (eg, diddleton and 
aungall, 1S52) require some knowledge of the physical 
characteristics of the sncwpack which may be difficult to 
obtain from either remote sensing or field measurements. 
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IV. Cevelo£{S6i)t cf a lecbnigae £cr tiodsling 
Inconiing Lccguave Baaiaticn in Mcuntainous Terrain 


One of the nest iitfottant parameters considered in the 
energy balance of a snovpaclc is incoming thermal radiation 
from the atmosphere and surreandiog terrain. This is a 
difficult parameter to seasure over an area, especially in 
mountainous terrain. It must be modeled very carefully to 
get meaningful results. There ace typically large diurnal 
variations and complicatlcns ace caused by cloud ccvecr 
forest cover, and sucreuxiding slopes. 

Longwave radiation is important in the energy budget of 
a sncwcovec because of the combined effects of the absorption 
cbaracteristics of scow in the thermal region and the 
attenuation chacacteristics of the atmosphere. Thermal 
radiation exchange with a snow cover is contingent upon three 
important considerations: 


1. Snow is a near perfect absernee of thermal 
radiation, while being a good reflector of solar 
radiation (US Srioy Corps of Engineers, 1956; 
Kondratyev, 1973). 


2. Clcud cover substantially affects the 
surface radiation balance by reducing incoising solar 
radiaticn and by incceasicginccmiag thermal 
radiation. This is especially important in high 
alpine areas where extensive cloud cover is likely 
during the snow season. 

3* Thermal radiaticn exchange cccurs during 
beth day and night. 

In order to realistically deal wirh the above 
consideratichs, the longwave radiation model should have the 
follcwing cbaiacteristics: 


1. It 
atmospheric 
day. 

2. It 
surrounding 
and thermal 

3. It 

cover and forest cover cn surface tbermai exchange. 

4. The mcdel should be driven by a limited set 
of discontinuous surface measurements, 

Work to date on the Icngwave radiation model has been 
limited to clear sky conditiens and does not include 
cortectiens fer forest cover effects. The assamption of a 
hoincqenecus snowcovereQ surface is made. The model has net 


should compute the variations in 
emissivity with elevation and time of 


should account fee the effects of 
slopes in terms of both horizon effects 
radiaticn of these slopes. 


should account for the effects of cloud 
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beeii;t^caplet( 'y tested at this tiisei though sufficient data 
haveleen collected (see Appendix 2)« The genecal structure 
of the model is presented in the fcllowing section. 

1. The study area is divided into a grid using 
the MCIC oigitial Terrain Tapes in a mannner 
described by Dozier and Outcalt (Appendix 1). 
MeasutemeDts of air teoperature and vapor pressure at 
various times and elevations are then interpolated to 
this grid using a lattice-tuning method and a 
bi-eubic spline interpolaticn. A tenperatuce and a 
vapor pressure is then stored fox each cell in the 
time/elevation matrix. This irfcrmaticu is then used 
to compute a teiiperature and vapor pressure for each 
grid point at each time step of the model. 

2. Thermal emissicn from the atmosphere is 
principally a functcn of the id ter vapor, carbon 
dioxide, and ozone ccntents (Paltridge and Platt, 

1S76) . Erutsaert (1975) develcped a tacnnique 
Mheceby atmosphereic emissivity can be determined 
from near surface vapcr pressure and air temperature 
measurements, eliminating the need for vapor pressure 
and temperature prefiles cr scuiidings. This 
technique is modified for the non-standard 
atmospheric conditens enccuntered in alpine areas. 
Ccrrecticns are made for both pressure variations, 
and temperature and vapor pressure lapse rates, which 
are determined from field measurements. Atmospheric 
emassivity is then cemputsa fen each grid point for 
every time step of the acdel. 

3. Qsing horizon information previously 
cemputed (as discusstd earlier) , a thermal view 
factor is computed for each grid point, using 
technigues developed by Reifscnxuer and Lull (1965), 
and Lee (1972), and modified by Dozier and Outcalt 
(Appendix 1) . 

4. Incoming longwave radiatou from the 
atmosphere and surtouding slopes is then calculated 
fob every grid point from the previously described 
information for a variety of time steps. The 
incoming longwave radiation matrix is then stored for 
later input to the snowmelt aicdai, and plotted as 
surface contour maps over the study area. 


Possible Errors and the heed for Further Jictk 

The weakriesses of the iiadel lie in its simplifying 
assunptions and in data measucGment errors, ’fne exclusion of 
cloud cover and forest covet parameters is respcnsible for 
most nedei errors. As reliable data becomes available on 
these parameters they will be included in the oicdel. Vapcr 
pressure data ^,crived from vet bulb tempetaruto is 


jsBisasss 


partic ilarly unreliable at cc around 0® c due to energy 
released during phase change. A sore reliable method for 
determining vapor pressure will be explored. Surface 
temperature of surrounding slopes/ used in Step 3 above ate 
generally poor estimates of reality. iJo account is made for 
slope or exposure differences and the actual measurement of 
the surface temperature is probably unreliable. Eetter 
techniques for measuring this parameter/ possibly using 
thermal imagery/ must he developed. In general a few sites 
where continuous measurements of all or the rnput parameters 
could be collected would greatly enhance the model's 
reliability. 
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V. Snow-Soil HsAt and Moisture Flow in the Southern Sierra 

The literature review on soil heat and thermal processes 
has teen essentially completed. Methods have been determined 
for estinating the uarious soil thermal and nydraulic 
parameters required for input ro soil heat-moisture transfer 
models. Varicus models have been surveyed and have Ipeen 
found to be both mathematically complex and therefore 
computationally expensive and partially inappropriate for the 
conditions under which the model will be used. 

Generally the literature has indicated that for the 
relatively coarse soils of the Southern Sierra a model of the 
^oil— snow heat and moisture exchange should have the; 
follcwing attributes; 

1. The model must be able to calculate heat 

flow by conduction with various surface ,bo=and^ty ) 

layer conditions Simula ting diurnal and, /"seasonal 

changes. 

2. The model shculd aGGount for downward 

percolation of meltwater# ceding the soil and 
thereby ameliorating upward heat t|an-sfer to an 
overlying sccwpack. | 

3. The model must adequately simulate the 
transfer of heat due to upward diffuii,^cn of water 
vapor# ccunter to the dewnward percolktion of 

meltwater. 3) 

'// 

■ // V- 

4. The model must also be useful in determining 
the formation of frost# the degree of ice segregation 
and the degree of ice saturatica with given initial 
cenditibrs, in order to evaluate the efficiency of 
meltwater runoff. 

5. The model must yield, as output# predicted . 
temperature and moisture content at various depths 
within the soil as well as an estimation of the total 
snow-soil heat and ncisture exchange. This allows 
testing and calihratcn using in situ measurements. 

6. The model must be compatible for use as a 
subroutine to the sncwmelt mcdel and its coupling to 
a hydrologic flow model. 

Many of the surveyed soil models require comracii input 
parameters. Keasurement and estimation of soil properties 
for input into soil heat-moisture transfer models ace 
presently underway. Samples obtained from tne thermistor/ 
soil mexsture block sites are used in these procedures. 
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Sail Eliysical 


Ercfetties 


Textural and structural analyses of six soil samples 
have been perfocaed using standard sieve-weight partitioning. 
These analyses determine the textural and structural 
character, both impcrtant for the estimation of thermal and | 
moisture properties.. The amcunts of clay and silt fracticns i 
in the various samples have been sufficiently low such that J 
particle size analysis of these fractions using differential 
settling methcds has act been used. In addition./to the 
standard textural and structucal analyses, mine rfa logical and 
predcainant particle shape icformation has been gathered 
using a standard binocular petrographic aicccscope with grain 
count grids.. 

■ The previously described measured soil pcopeties ate 
used in the calculation cf soil pcccsity, bulk density, 
optimum moisture content and ether 'parameters. These are 
used in the estimation of soil thetmai and maisture 
properties. 


The soil thermal conduct 
estimated by two methods. Te 
compesitien, and particle-sha 
samples will be used as input 
emplayed by Philip and DeVrie 
properties. These estimates 
experimen ta 1 va lues cf the r la 
for similar soils as determin 
temperature dependence cf the 
capacity due to the presence 
separate of soil (Penuer, 197 
calculations with our sell sa 
fines present. However these 
change wit b liquid water and 
Hcacberts, 1973) . 


ivity and heat capacity will be 
xture, minetdlcgical 
pe distribution of the various 
to the material f racten ocdel 
3 ( 1 557) to estimate these 
will be campared with the 
1 cendcctivity and heat capacity 
ed by Kersten (1949) . The 
rmal conductivity and heat 
cf unfrozen water in the fine 
0) is neglected in the 
m pie s due to the lo w q a an tit y of 
parameters do effectively 
ice content ( Ni xon an d 


The hydraulic properties of the soil samples will be 
determined using textural and particle shape information with 
the appropriate moisture characteristic and desGrpticn 
curves. Outcalt (1976) specified icg-log fit parameters to 
abstract the shape cf both the desorption and the unsaturated 
hydraulic conductivity curves. Soil porosity was used with 
these parameters to establish the relaticnships between the 
hydraulic head, the unsaturated hydcauiic conductivity and 
the water volume content. These input values can also be 
used to describe the therffal flo# as it is; coupled 

with the moisture transfer. Green and Corey (1971) have 
reviewed varicus methods of calculdting unsatarated hydraulic 
conductivity using sell pore-size dlstcibutuion solely. They 
concluded that all of the reviewed Kethedsadeguately 
predicted experiientaly measured values for a range of soil 
types. To obtain estimatiens fob the soiia of the Sierra^ 
Nevada , more than one methed will be used to calculate each 


of the 30x1 hydraulic properties. 
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Fiela 2l€asur€[0£ats 

Durinig the smaner and fall, 1977, field temperature and 
moisture measure men t sites here established in the test area. 
The locations of the test sites were chosen on the basis cf 
ease of access, Iccatcn of the Cccperative Snow Survey snow 
courses, and the representativeness of the locaticn. 

Seven thermistors and six soil moisture blocks were 
installed at each site at various depths. These were 
sel'A^cted to correspond to textural or composition changes 
obsffrved during the installation. In some sites a total 
depth-range sample was prevented due to difficulty in digging 
holes. At twc sites, twc soil moisture blocks were assigned 
to each depth with separate holes. At the ether sites all 
six cf the mcisture blocks were installed at various depths. 
There were three thermistors placed within the soil at each 
site at depths corresponding to these of some of the moisture 
biocke. During the installation soil samples were obtained 
from the appropriate depths, and a preliminary field 
desetiptioa was done. 

Since the test sites were established eacn site has been 
visited during the scow season to obtain in situ measurements 
of soil and snow temperature and soil moisture content (see 
Apperdix 4) . 
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Soil Moisture Elcck Accuracy 

The soil moisture blocks measure resistance across 
stainless steel grids implanted in fine- textured 
plaster-of’-paris (gypsumy blocks. The boisture measured is 
transferred from the surrounding soil into the block by the 
capillary suction of the block material. Tae resistance 

corresponds to the moisture content in 
teriiis cf percentage cf available water for anfroxen 
conditions. Although it is reascnable to assume that the 
block measureiients ace accurate due to the textural contrast. 
between the block and the surrcanblng soil for unfrozen 
conditions, it has been shown that the nlock resistance 
increases greatly when freezing occurs (Sartz 1967,1970). 

Also it is expected that the soil water system would freeze 
and thaw sooner than the block water system due to their 
respective textures. In order to calibrate the block 
ffleasurementg airing freezing, other moisture seosets have 
been obtained with diffecent ccristruction. Tuese sensors 
operate on the same principle as the blocks, DUt cersist cf a 
sandwich of stainless steel and fiberglass. Toe porosity and 
perineibility cf these Sensors is siiriTar to tnat of the soil. 
These sensers contain built-in tkermistors and may be 
oriented such that they are most sansitive to vertical 
moisture moveinent. They will previde contrei data fee 
comparative use with the soil moisture block data. 
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VI. Gecgcaphic Oesciiptioo or the Fergana Area 


Intrcduct icn 

This study pcovides a tasic geographic description of 
the Fergana Valley in the South Central OSSft, with emphasis 
on the Soviet study area, part of the joint US/OSSB 
cooperative research program, which includes rhe basins of 
the Gavasay, Sumsar# Kassansaj and Eadsnaata Hivers-^ It 
contains information provided by the^^ Soviet Unicn on the 
study area, as well as infcrmaticn contained in earlier 
progress reports. 

Soviet Central Asia is part c£ a vast synclinal 
depression between the Bussian and the central Siberian 
Platforms. Bolling plains, typical of the region, are made 
distinctive by the effects of a desert climate, while the 
southern rim is accented by the Pamirs, the highest 
mountains in the QSSH. The Fergana Valley (Ferganskaya 
Oolina) , a large intermcntane sy nclinorium, lies in the 
central eastern part of Scviet Ceritral Asia (Figure 1) , just 
north of the Pamirs. The Kucamin and Chatkai mountains 
(noth part of the Tien Shan Bangs) bound it to the north, 
the Alai and Turkistan to the south, and the Fergana 
Mountains to the east (Figure 2). The floor of the valley, 
sloping frcm east to west, is separated from the foothills 
of the surrounding mountains by a series of weli-developed 
alluvial fans. Flowing westward along the northern edge of 
the valley is the Syr Darya. South of the Syr Darya is an 
area of relatively level land covered by sand dunes, salt 
marshes, and salt lakes. A large alluvial fan at the 
southwestern end of the valley was formed by sediments 
deposited by the swif t-flcwing Sckh Biver, which disgorges 
into the valley from the mctntains in the south. A zone of 
springs and wells occurs around the perimeter of the 
alluvial fan between the 430 and 400 m contours (Plummer, et 
al, 1S71). Cn the west is the narrow (8 km wide) Leninabad 
Pass, through which the Syr Darya flows, joining the Fergana 
Valley to the Golcdnaya (Hungry) Steppe. The Soviet study 
area, in the north-central pctticn of the valley, is bounded 
by the Kuramin Eange to the west, and by the Chatkai Bangs 
to the north and east. Rivers, w supply the Syr Darya 

drain the southern slopes of these mountains (Figure 2). 

^ Fergana Valley is at roughly 40^ to 42^ N 

latitude, approximately the latitude of Denver, Cclcradc, 

USA iFigure 3)i with a longitude cf 70° to 73°. The 
valley is in the font ox an ellipse, 300 km long; and 150 km 
wide, with an area of acpccximately 22,000 sg km. Thsi 
surrcunding acuntairs reach an elevation of nearly 6000 ra in 
the south and 4500 m in the ncrth while the surface of the 
valley rises frcs 300 tc 4CG m in t west (Leninabad is at 
325 II) to 500 m at Andi 2 har. in the eastern end, and 900 to 
1000 m in the southeast (Csh, in the foothills, is at 990 m) 
(Suslcv, 196 1). JiaBangan, in the north, is at 440 m, and 


Ferqana, ic rie south, is at 580 o. 'The surface cf the 
valley is fairly level, the central part of which is desert. 
Salt playas, for the most part dry, but flooded in some 
years by high river water, are cbaracterisric of this 
central region, as are salt-aarsh swaaps, such as the 
Dam-Kul, which are in the prccess of drying out (Suslov, 
1961). 
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The Fergana Valley lies mainly in the Uzbek SSa and 
partly in the Tadzhik and Kirgiz SSE's (Figure 4). (The 
intricate coundary delineation in the area stems from an 
effort to include the acuctain watershed areas with the 
parts of the basin to which they supply water, and 
considerations cf different nationality groups)* Although 
Uzbekistaa contains only cne-tbird of the total area in 
Central Asia, the Fergana Valley provides it two-thirds of 
the population, three-fourths of the industrial output, and 
over ene-oalf cf the arable land. 


The Fergana Basin (Figure 5) has the highest population 
density and the highest ccncentraticu of agriculture and 
industry in Uzbekistah (Gregory, 1968). The foothill ring 
is the most densely populated; here there is a continuous 
belt cf cultivated, irrigated land which receives water from 
numerous 'mountain st teams. This large population 
concentration in a small area contrasts with the large and 
virtually uiiinhaiited desert and semi-desert areas cf most 
of Soviet Central Asia (Allwcrth, 1967) . Most of the 
Fergana Valley is characterized ty a rather dense rural 
popolaticn of mere than 50 persons per sg jua, rural 
population making up 60^ or mere of the total population of 
the valley (CIA, 1974) . Howeyer, the western end cf the 
valley, in the Kckand region, has a rufel population of only 
1 to 10 persons per sg km. The largest cities in the valley 
are Andizhan (172,000) , Haaangan ( 164,000) , Kokand 
(133,000) , Fergana (11 1 ,000) , and Rargelan (89,000) , 

Kuvasai and Leninsk haye pcpulaticns of 30,000 to 100,000. 
Several cities alone the scuthera edge of the valley have 
populations of 10,000 to 30,000, with many smaller cities of 
les$ than 10, 000 (Eurikin, 1967) . The guantity of water 
that is available partly esplainsth^ density of villages 
located on the lower elevations alluvial fans and on 

the numerous distributarie,'-- cf rivers (Plummer, et 

al., 1S71). , \ / V, . 


The vegetation of the Fergana Valley (Figure 6) forms a 
series of belts, or tings, scrrcuiiding the central Fergana 
desert and extending into the foothillk' and mountains 
surrounding the valley. The nature of the plant cover of 
this territory varies widely in relation to cfae local 
elevation, slope expesure, and cdinfall distribatuicn . 
Several main vegetative belts can be distinguished: 1) 

desert, 2) steppe, 3) itcurtain-meadow, alpine and subaipine, 
4) forests, thin-ferests, and scrub, and 5) rocKy ridges, 
mounds , moraines, and firns with sparse vegetaticn . 



20 




m 




m 


« 



The vegetation of the fccthill desert is oainlj 
represented by ephenierals, scnetiiies in community with 
ucrinwccd and Russian thistle. Wcrmuccu is round cn the 
alluvial fans and into the foothills. A tyt»ical feature of 
ephemeral vegetation is the fact that it develops ctly in 
the vet autumn feriod and estecialiy in the spring period, 
when the grass stand is continuous, even though shcrt-lived 
(30-40 days), Sith the cessation cf the spring rains, the 
above-ground part of the ephemers rapidly dries up and the > 
locale changes into a sun-scorched gcay-yellov desert vith 
individual bushes of ucrnijccd or Eussian thistle, carpeting 
only 30-40 percent of the ground. 

Steppe areas are chdracterized mainly by 
drought- resistant vegetation represented by various cereals 
and some spring-autumn ephemecs, 

aountain meadovs are chiefly ccmpcsed of ” isescphy tes” , 
which require year-round precipitatiou for their grcuth. 
These mehdous are found at slevatlcns of 1500-1300 m. The 
grass stands typically grew to 1 m and 80-50 percent of the 
soil is cowered by vegetation. s 

Subalpine vegetation is found at eievatiens of 
3000-3500 m. Grass stands are dense, covering 85-$C 
percent of the soil, and it is 35-50 cm high. 

Alpine meadows have the highest location (3200-4000 m) 
and ace represented by motley grass-cereal ana sedge ■ 
meadews. The height cf grass stands here are 25-35 cm and 
the soil is 50-80 percent covered by the vegetation. r 

Forests, as commonly understood, are scattered in small 
sections. Spruce, spruce-fir, juniper, and nut-bearing 
forests are most common. Among the trees are found an 
understory with infieguent ciusters of spirea, dog rose, and 
honeysuckle. o 

Reeky sectiens of ridges and slopes, mounds, and Recent 
moraines at the highest altitudes have sparse plant cover. 
This zone has roughly the same appearance over the entire 
territory described. 

Table 1 lists the distributicn of vegetative types, 
by percent area, in the test site tasins. 
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Figure tJ"- Comparative area and latitude of the 
US and USSR (from CIA» 1974). 
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Table 1. Rock complexes, soil cover and vegetation of test area baslns=^(by percent area) 




Clioate 


The Fergaaa Valley lies in a semi-arid regionywith Ccld 
winters and dry, hot summers. The Koppen-Gieyer climate 
classificaticn fct the region falls in the boundar!^ zone 
between a dry, cold, desert (EWJs) and a dry, cold steppe 
(BSk) (Strahler, 1974) . The Koppen System classified the 
desert regions of Sewada-Ctah in the western United /States 
and the Atacama Cesert in Chile as having similar climates. 
This ccmpariscn is deceptive. The climate of the fergaaa 
must not be viewed cn a macrcscale, nut a mere mesc cr micro 
level to fully understand the impact of climate on the 
valley’s agricultural and industrial productivity. 


The two most important physicgraphic factors in the 
climate of the Fergana Valley are its continental location 
and the surrounding mountain ranees. The valley’s 
qeaqraphic location basically remeves any marine effect and 
produces relatively extreme teiuperatures in all seasons. 

The surrounding mcuctains modify this somewhat by providing 
a rain shadow effect in some seasons and intense crcgtaphic 
conditions in ethers. In additicn, the surrounding 
mountains give the Fergana yalley a more mild climate than 
it might have, although snow and frost do occur in the 
valley preper in December and January, when the temperature 
may fall as lew as -15®C. At Tashkent, a first order 
weather station about 80 miles nerthwest of Leninatad and 
outside the Fergana Valley itself, temperature wariatiens 
are much more extreme. This may be seen in the mean monthly 
temperatures, averaged ever the period 1951 to 1960, shewn 
in Table 2. Figrte 7 shows climcgraphs of Fergana, 

Namangan, and Osh (Walter and leith, 1967). 

Soviet Central Asia is basically under the influence cf 
travelling westerly air masses, with the cycio- and 
anticyclogenesis intrinsic tc it. Also quite important are 
cyclcnes, moist westerly air masses, and cold northerly air 
masses, lowering ait temperature and causing precipitation. 
The ' relative abundance of spring and winter precipitation is 
caused by cyclcnic storms from the Mediterranean and Elack 
Seas. In summer, the air is too dry for precipitation in 
spite cf or og ra phic cu t r en t s . T h e mo unt ains to the so ut h 
provide a barrier prcducing a rain shadow effect from the 
mcnsccn of Southern Asia (Sc slov, 196 1). 

Owing to the divetsity cf slope exposure with respect 
to the moisture-bearing air aasses, precipitaticn varies 
thrcuqhout the regicn. The ttean-annual precipitaticn totals 
fluctuate from 80-100 to 4C0-750 mm. (Resources of 
Surface Waters of the fJSSS, 1969) 

The sma 1 lest annual precipitation totals are otserved 
in the Gavasay Diver basin, cccupying^^ most westerly 

position in the study area. As cce moves eastward, 
precipxtar ion increases and reaches maximum values in the 
Padsbaata River basin. On the soi'theast slopes of the 




Chatkal range little precipitaticn fall:.. In the valleys of 
the Gavasay and Kassansay tivecs, at 1000-1800 B,tne total 
precipitation each year dees not exceed 300 om (aesources of 
Surface Waters of the USSB, 1969). Heavy precipitation 
falls on tie east side of the Chatkal range. In the 
Padstaata River basin, at 1500 a elevation, the annual 
average is 700 ma; and near the crests of the acuotains, the 
average is 1000-1200 aa. Figure 8 shows the distribution of 
annual precipitation totals with elevation in this tegion* 

In the foothill and low-mountain zones, most precipitation 
occurs in winter and spring (Figure 9); summer and autu|mn 
are dry. Table 3 gives mean aenthly precipitation and 
water balance data fer Tashkent, Fergana, Andizhan and 
Namangan. 

Cempared with the focthills, precipitation totals are 
more evenly distributed throughout the year in the aountain 
regiens. At about 2000 m there is little difference between 
sumaer and winter prccipitaticn totals. f\ 

In most cases, in the acuntain region, the number of 
days with pracipitation goes up with a^ concommitent increase 
in precipitation totals: the aeah annual number of days 

with precipitat ion in this region is 100 to 150 days. 

Some years, precipitaticn in the form of rain showers 
occur five to seven times a year in the basins of the 
Padshaata and Kassansay rivers. The showers may last 20-25 
hours. Maximum sbeset intensity varies f too 0. 10 to 1. 60 
ma/min or higher. In June 1957, the aaximua shewer 
intensity in the Kassansay liver basin teached 3.1 mni/min 
(Besources of Surface Waters of the USSB , 1969) . 

Owing to the gecgraphic positen, this region gets a 
relatively large amount of sclar energy, especially in 
sumaer. The influx of solar radiation in the foothill part 
of the basin in a year’s tine is about 140 kcal/ca^ 

(Hanebeok on Clinate in the USS8, 1973). With increased 
elevation, owing to a reduction in the density and dust 
content of the air, the total influx of radiaticn rises to 
150-155 kcal/cm2 (Handbeok en Climate in the USSl, 1973). The 
largest mean-monthly levels cf total radiation are about 20 
kcal/cvz and cccur in June and July. The smallest 
radiation influx is observed everywhere in December; it is 
4-7 kcal/cm^ (Handbeok on Climate in the OSSH, 1973). 

Annual levels c£ net radiaticn change bver the 
territory from 46 tc 50 kcal/cm^ (Spra voennik po klimatu 
SSSR, 1973), The laEgest total fficnthly levels (8-11 
kcal/cm2) cccur in June and July (Spravochnikpo klimatu 
SSSB) Net radiaticn is Icwest in Cecember and January, 
when it drops to 0.3-0. 5 kcal/cm^. 

The relationship of evaporation with local elevation is 
manifested guite clearly. The ameunt of evaporatiGn depends 
primarily cn climatic facters; mainly precipitation and 
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relative huibiaity. Ihe latter decreases abcuttly vith 
elevation in the ncuntains. a graphic presentation of th 
decrease in vapct pressure eith elevation is given by Fig 
10. Because precipitation increases with elevation (to a 
certain point) evaporation also initially increases with 
elevation; at higher elevations, bo hevec, it decreases 
(Shul ' ts, 1S65) . In the valley proper, in sunnier, clear 
skies and dry sukhovey winds from the southeast can cause 
evapcraticn which may exceed annual precipitaticn by many 
times. 


e 

ure 


In the mcuntaincus part of the basins, air temperature 
decreases with elevation. Fcr example, at eievaticns to 
1200 m the mean daily air temperature in January is frca -2 i 
to -50; at 1200-2000 ra, -3 to -7^; and in the 200C-300C m 
zone , the mean daily air tempstat ure drops to -d to -9° 
(Besources of Surface Waters of the USSfi, 1S6S) . July and 
August are the hottest months. The highest mean daily 
tern peratures ars observed in the foothills (2b-28°) . At 
the elevations of 2000-3000 m it 'does not exceed 1C-20®. 

In the high-mountain zone night frosts are common in summer. 

In winter, on the slopes of the Chatkal range the 
adiabatic lapse rate ranges f rom 0^27 to 1. , and in 
summer gradients range from 0.11, to 0.73° (Besources of 
Surface Waters of the DSSH, 1S69| , 


0 


0 . 

K 





Table .2.- Mean Monthly Temperatures 
for Tashkent, 1951-1960 



°C 

°F 

January , 

.9 

33.6 

February ^ 

2.8 

37 

March 

6.3 

43.3 

April 

14.3 

57.7 

May 

19.5 . 

67.1 

June 

24.5 

76.1 

July 

26.9 

80.4 

August 

24. 8i 

76.6 

September 

19.5 

76.1 

October 

12.4 

54.3 

November 

5.3 

41.5 

December 


34.2 

Mean 

13.2 

55.8 


(from World Weather Records, 1967) 
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Figure 7. CHmographs of Fergana, Namangan, and 
Osh (Osch) (from Fetisova, 1965). 
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ST 
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66 

72 

61 

44 

27 

15 

9 

7 

6 

14 
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AE 

0 

0 

20 

30 

38 

28 

18 

9 

4 

15 

9 

0 
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0 

0 

0 

0 

37 

78 

136 
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95 
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19 
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9 

8 

12 
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0 

1 

22 

45 

43 

31 

19 

9 

6 

20 

13 

1 
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0 

0 

0 

0 

54 

77 

.128 

ISO 
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82 
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0 
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2 

21 

61 

115 

160 

180 

150 

86 

42 

14 

2 

833 
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45 

37 

63 

51 

29 

12 

3 

1 

5 

26 

34 

42 
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ST 1 

127 

162 

204 

197 

148 

90 

50 

30 

23 

22 

42 

82 


AE 

0 

2 

21 

58 

78 

70 

43 

21 

12 

27 

14 

2 

348 

0 

0 

0 

0 

3 

37 

90 

137 

129 

74 

15 

0 

0 

485 

S 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 . 

0 

0 

0 


PE Potential evapotranspiration AE Actual evapotranspi ration 

P Precipitation D Water deficit 

ST Soil moisture storage S Water surplus 

(from Average Climatic Water Balance Data of the Continents, 1963) 
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Sells 


Soils of th<3 Ferjana Valley 
eater-physical ptcpeities. Xhe s 
recede uinate in the region, in con 
'^nd parent material, have produce 
Sierezem* Sierezem soils develop 
uindhlown deposits (loess), along 
alluvial fans where rivers: emerge 
soils cover only 2% of the USSH, 
valley (Figure 11) . Small areas 
and lithoscls can be found inters 


have widely varying 
emi-atid conditions that 
junction witn topography 
da gray soil knewe as 
under arid conditions on 
piedmont areac, and on 
from mountains.' Sierozeo 
tut dominate the Fergana 
of soils such as Regcscls 
ittently in the valley. 


In the foothills to 160C m elevation, siecozems develop 
beneath a grassy cover with ephemecs (Figure 11). Ihe 
sietezems have a scaly"^laatinat6d structure in the surface 
layer, changing into a siightiy cloddy layer. Below lies a 
weakly compacted herizet with segregation of grains and 
concretions of catbenates, underlain by a layer rich in 
carbonates and gypsua. Biogenic constitution 
(f oraminatedness and vesicularity ) is distinctly evident in 
the sierozems of different kinds. Ihe carbonate profile is 
well-defined, even though the entire soil is carbonate in 
character. The content of ncn-silicaceous salts is high, 
from 5 to 20 percent or more of the total soil mass. The 
sierozems differ in content of elements with waten-resistant 
micrcstructure, that is, structural entities larger than 1.0 
mm. But on the other hand, the water-resistant 
micrcstructure is well-defined; in fact, it is determined 
by the fairly favorable water-air properties - good water 
permeability and high, mainly capillary, porosity. The 
humus content in sietozeins fluctuates from 1 to 4 percent. 
Ihe specific weight of siecozems is 2.5 - 2.8; the bulk 
weight is 1.0 - 1.6; and t he porosity is 50-60 percent. The 
total water permeability of the sierozems is 0.5 - 2.1 
aa/min. (Resoiirces of- Surface 'Waters of the USSH, 1969) 

Characteristic of the sierozems is a ncnleaching type 
of water regime. The soil layers begin to ne wetted in 
autumn or winter:. Total wetting depth varies from 1 to 
2m. Sierozems in the western Tien-Shan in a 2 m layer have 
a water reserve cf 382 mm (Karpov, 1960). At the end of 
Hatch or the ceginning of April vegetation starts rapidly 
drying up. Dessication is usually completed by May (Bode, 
1963) . :■ ^ : ■ : 

Osing classif icatioii schemes develGped in the C£, 
SierGzem soils can be classified endec the order of 
Aridisols or Xercsols. The classif icaxicn Aridiscl is from 
the Cempretensi ve Soil Survey System and Xerosol from the 
FAC-CKISCO soil map and ’•Soil Systetr of the Wcrld'’ {3rady, 
1574) . : ■ . ■ : 

Arifliscls are mineral sciis fcand in arid climates. 
Except where there is greurd water or irrigation, as is the 
case in much c£ the Fergana Valley, ail of the soil horizons 
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ace dry moce ckan 6 aontbs of the year. Conseguently , they 
have not been subject to intensive leachiny. They have an 
ochric epi^edion genecally light in Color, loW in ocganic 
fflattec, and may be hard and sassive when dry. They aay have 
a hcrizon cf accuaulaticn c£ calciua cacbonate (calcic) , 
often with an underlying layer of gypsum , 3 or even acre 
solutle salts (Brady , 197*1) - The Sierozea soils of the 
Fergana Valley come under the subcrdet Orthids. This 
suborder is broken down into t»o groups. The first group is 
coopcsed of soils of gently cx incderately sloping terrain, 
used mostly as range land with soae irrigated crcps. Areas 
in the United States that have 
northeast Arizona, central and 
southeast eolcradc, and scuth-ceotrai Idaho. The second 
group includes soils of gently sloping to steep terrain. An 
area in the US that has a similar soil is norxh-central 
Nevada. Orthids are differentiated from other subocdecs of 
Aridcscls by the absence of a horizon of clay accuaulation. 


3016 irrigated ctcps. 
similar soils include 
ncrxbeast New Bexicc, 


Xeroscl, another term used for the classification of 
d'eseit or semi-desert soils is broken into four suborders, 
only three of which apply to soil types found in the Fergana 
Valley. They ate haplid (typical xecosols or semi-desert* 
soils) , calcic (semi-desert soils with h distinct lime 
accuttuiaton hcrizcn) , and gypsic (semi-desert soils with a 
distinct gypsum accuiiiulaticn |hcrizcn) . 

Sierozems are typically rich in calcium and form gcod 
agricultural soils iVhen irrigated and fertilized. aithout 
itriqaticn, Sierozems are not suitable for growing crops. 
Areas that are not irrigated are used for sheep and cattle 
grazing,. but prodcctiwity per unit area is low. 

Brown and chestnut soils are found at hiyher 
elevations. These soils develop under sheep's fescue and 
gjuack qrass-mctley grass-brush steppes, and can be found to 
an elevation cf 3200 m. Loose loess and loessiike rocks on 
deluvial and eluvial deposits are the soil-forming rocks of 
this group. In mechanical composition they are 
fine-grained, rubbly, and rubbly-gravely. Humus content 
varies from 4 to 20 percent. The structure of the huiius 
horizons is highly vater-resistact and water-permeable, 
nutlike, or cloddy-granular, intergrading into ccarse-cloddy 
with depth. The carbonate content in the upper part of the 
profile is 1 to 4 percent, and in xhe lower profile - to 10 
percent. Physical properties of these soils may be 
described as follows: specific weight 2.5 - 2 . 3; bulk 

weight 0.6 - 1.1; and porosity 56 - 66 percent. The water 
permeability is 1.8 - 2.1 am/min (Fesources of Surface 
Waters of the USSR, 196S) . 

light-br cTni meadow -steppe so 
secticns of the slopes of the ran 
reaching the perpetual sncws. Th 
from 3 to 20 percent humus , are w 
horizons, and are rubbly. Soils 


ils cccupy the highest 
ges and their watersheds, 
e soils are th in , contain 
eakly differentiated by 
underlie suoalpine a nd 


alpis6 compact oeadcu-stoppe ij£as£ stands on deluvial and 
eluvial Paleozoic aod Bccaine deposits* Sodding oz the 
uppex horizon and good gcannlac structure is typical. The 
total porosity is 45 - 61 percent; porosity in the huous 
horizon is 57 - 60 percent. 

Precipitaticn exceeds evaporation in the zcne ci 
Deadcw»steppe soils; this ensures a leaching type c£ water 
regiae. Mcst of the zone is characterized by near-surface 
bedding of scil-grcunduater ; its surface rises into the soil 
profile during snow lelt, during the spring-aut uan tains, 
and scaetiaes even in heavy suaaer rains. 

In the Padshaata Bivet basin the light-brown 
aeadcw-steppe hig'h-acuntain soils (coarsely skeletal) ate 
developed; they fora aaidst cctcroppings of bedrock and in 
the region of glaciers and firns beneath sparse aeadcu 
vegetation, and are to be found in small isolated patches. 
The soil body is continuously wetted. 

E-rincipal soils in the Soviet study area are listed in 
Table 1. 

A major ptoblGa in the agricultural areas of the test 
site is soil salinity. In areas cf prolonged cultivation, 
such as the Fergana, excessive salts may accamalate. This 
is especially true of the central pcction of the valley. 
These soils reguiie careful drainage and salinty ccntrol. 

The soils near the perifteter cf the valley, including the 
foothills and lever slopes of the meuntaias, are reasonably 
unifera in physical and cheaical composition and do not 
contain harmful salts (Allvocth, 1967, Burikin, 1 967). 

In the central and western pacts of the valley there 
are sands and alkaline aarshes. At the bases of the 
mountain ranges, broad fans of preluviai and deluvial gravel 
deposits occur (Suslcv, 1961). i Soils in these areas are 
not cf much agricultural importance. They fall under the 
great soil groups of Eegcscls and Lithosols. Hegoscls are 
weakly developed soils cn Iccse, non-alluvial deposits 
(loess, sand) . The three suborders that apply are Butcic 
Begcscls,-%alcaric Begosois, and Bystric Begosols* 

Lithcscls are ncn-acable scils, not even fit for grazing 
under most conditions. They ace primitive, cubbly soils no 
more than 30 cm deep on bed rock (Ecady, 1974). 

Suslov, writing on the central desert pertien cf the 
valley states; "Intensive processes of weachering of 
Tertiary conglomerate' and sandstene, and also of recent 
sandy alluvial deposits alcng the Syr Darya, lead to the 
creation of ccnsiderdble areas of shifting sand. Cnly along 
the margin of the valley is the sand somewhat secured by 
vegetation. The partially mebile sand of tae center 
accumulates in recent barchans, with heights up to 50 feet, 
which move east, covering fields and settlements." 









fcFiting about the gravel detcsits, Suslov makes note 
that the gradual ercsional "disintegration of these 
formations creates, at first, individual table mountains, 
and, later , a highly broken relief. The rivers have large 
d dry deltas, as for example, the huge alluvial fan of the 

Sokh River, which irrigates Kokand Casis, and the delta of 
the Isfaytam River." 
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Completion of cap ti oh to Figure 11 on preceding page 

% 

01. gray-brown desert, takyr, and sandy soils of the plains 

2. zone of sierozems, mountains and foothills 

3. zone of dark brown and chestnut soils 

k. zone of 1 ight -brown soils of high-mountain steppes 

5. high-mountain and desert-steppe takyrlike and syrt soils 

6. strongly salinated soils on salt-bearing Tertai ry deposits 

7. light-brown, meadow-steppe, high-mountain soils, eroded, coarsely skeletal 
amidst outcroppings of bedrock, firns, and glaciers 

8. boundary Of region 
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Geology 

The Fergana Valley lies on the north aestecn edge of 
the Edfflir Range ^ thel highest acuntain range in the CSSR. 

The Eanirs «ere upli(ft6d fchec the Indian plate collided with 
the Asian plate during- the Cenozoic (Dietz and Holden# 

1970). Tectonic activity frcis plate moveinents during the 
Quaternary created the Fergana Valley and the ocuntain 
ranees which surreund it (Figure 12) . The Fergana is a drop 
blocX valley surrounded by rapidly uplifting mountain ranges 
(Eelcusov and Gzovskii# 1958). The area is descrihed by the 
Hinistry of Geolcgy# USSR, as a region of perigeosynclinal 
conplexes which are part of the tectonic regime of the 
Scythc-Turaiiian structures (Ministry or Geology# 1964). 
Pigure 13 is a stylized map of the Fergana area shewing the 
general teGtcEic structuce and seisnic activity. 

The Fergana gccsynclinal regicn can be divided into 
mountain and foothill sectiens. The study area includes 
diverse topography - from alpine zenes through regiens of 
oediuni- and Icw-elevaticn# tc footnills and flatland. The 
average elevatidns of the river watersheds in the study area 
vary from 2300-2800 b. Regions of mediutt-elevaticn relief 
are dissected by deep and narrow valleys# while regions of 
low-elevation contain numerous valleys# gorges# ravines and 
kettlehcles. 

The geologic structure of the study area and its 
surtcundinqs are associated somewhat with the terrain. 

There is a ccnsistent shift from the periphery to the crests 
of the ranges of geologic hlccks which differ in age# 
compesitien, degree of dislccaticn and aydrologic 
cha tacteristics. Oldet# ccnsclicated igneous and 
metamorphic recks ate found at higher elevations while less 
consclidated Cuaternary deposits are found at the foot cf 
slopes, ale, ng the fleers cf, valleys# and in karst tasins. 

The Kuramin and the southwest porticn of the Chatkal. ranges, 
are eomposed chiefly cf recent iDteusive igneous rocks and 
older (Paleozoic) metamcrphic recks. lutrustive and 
extrusive igneous cocks contain fissures more than 50 m deep 
and hence for B geed water xesecvcirs. Clayey and sandy 
slates# also widespread# are nearly impermeanie . Fissures 
contained within them ate usually narrow <ind shallow in 
depth (to 30 m) # so their water stcrage capacity is alSG 
small and tew springs issue from these rocks. Pure 
limestones are r:arely found in the area. In i-laces# 
deluvial fcrniaticns cverlie the broken terrain of the 
mountain ranges. The foothills and valley floor are 
compcsed of Tertiary sediments and v^Jaxternary alluvium from 
the surrounding mountains, lithclogic ccaplexes in the 
watersheds are shewn in Tab le 1 ^ 

The Fergana Valley is selsmicaliy active. A major 
fault runs al eng the i north side Cf the vaiiey at the base cf 
the Kuramin# Chatkal and Fecgana Kcuntain Ranges. Faulting 
is evident on all Sides of the valley fram geophysical 
evidence and gectbacadl activity (Ministry af Geology# 

1966E ; Earanov, 19S9) . 




Figure 12. Fergana Valley and surrounding mountain 
ranges. 
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Figure 13. General tectonic structure and saismic 
activity of the Fergana Valley region (after Belousov 
and Gzovskii, 1958). 


Hydrogeolcgy 


Grcundwater conditiocs in the Fergana Valley are the 
result of varying aguif er-sysrem characteristics such as 
aqhifer materials/ altitude/ Iccaticnal relation to 
irrigation outflow points/ and location in relation to other 
water-bearing or water-confining rocks. The hydrcgaology of 
the study area is determined largely by the structural 
geology. Due to the varying water-tearing nature of the 
different rocks and unconsolidated deposits and their 
subsurface spatial variation/ the groundwater 
characteristics vary spatially even within a gi-ven water 
body (figure 14). 


Infiltration in the mountainous areas supplies 
groundwater to the Fergana basin. Inere is however little 
surface water infiltration at the higher elevations in the 
zone of igneous and iretamcrpbic rocks. Infiltration mainly 
occurs from channels which flew through the less 
consolidated late Tertiary and Quaternary fluvial/ alluvial 
and glacial deposits. Host infiltration replenishes 
groundwater; some re-emerges in depressions both at 
springs and riverS/ and some supplies flows in the alluvial 
fans of the streams. Groundwater in the study area is 
supplied primarily by infiltration from surface streams. 
Carbonate bedrock provides a supplementary source of 
groundwater. Atmospheric precipitation typically plays a 
small part in supplying the groundwater/ since the water 
table lies at considerable depth. The denudation of the 
mountain ranges and the large local slopes provide good 
drainage and favorable conditions for intensive water 
turnover, a.ecause of the slopes in the Zaadyr and Hezhadyr 
depressions/ groups of greunewater flows/ supplied by 
infiltratcn from nuiietocs surface streams/ are formed in 
this part of the study area. 


I The aquifer sys 
include three categc 
uncenfined and semi- 
water body. A shall 
that occurs at shall 
perched water body) 
groundwater bodies, 
supported by a relat 
sediBents, silts or 
younger alluvium and 
where water bearing 
lac us trine an d m ars h 
result in salinity p 
of the Fergana Easin 
d ev el oping n c ti z c n s ) 
alluvial terraces (S 


tern of the a 
ties: 1) sfa 

confined wat 
cw water bod 
cw depth in 
that is usua 
Shallow wat 
ively imperm 
cla y. These 
probably in 
deposits are 
sediments, 
rob lets i and 
in areas of 
/ and are co 
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rea may be broken down to 
allow water body/ 2) 
ec body and 3) confined 
i is defined as groundwater 
shallow wells (this may be a 
lly semi-discrete from other 
er bodies are comracnly 
eable layer of fine 
bodies occur only in 
the flood basin deposits 
i n terla yered with 
These water bodies often 
mainly occur in the bottom 
marsh and solonchak (or 
mine a in the Adyr area of 


Unconfined water bodies may be identified frero w^^^ 
data as having little cc no difference in heau with-^^v 
depth. These usually cccur in unccnsciidated deposit s that 






ace not intecteoded vitb fine (jcained lacustclne cc aacsb 
deposits. Tertiary ccnglocerate and saadstone overlain in 
the foothills by broad fans of prcluviai and deluvial gravel 
deposits such as in the desert pastures on the margin of the 
Fergena Basiu arc conducive to such environments, Ihe 
unccnfined qrcundtiatec bodies found in the IZaadyr and 
Bezhadyr depressions of the study area have no pressure head 
and lie in boulder** gravel deposits at considerable depth. 

Semi- confined and confined conditions occur vhere the 
water body is bounded at the base and/or surface by 
impermeable layers of soil, cegolith, or rocJc* These water 
bodies tend to occur in Icwer valley areas adjacent to 
foothill terraces and in the central portion where alluvial 
materials, loess, and vacyin g-facie fluvial material are 
interbedded and interf ingecing. In the Soviet study area, 
they occur along the periphery of the Zaadyc and Mezhadyr 
depcessicns where water-bearing bccizons are enclosed 
between relatively impermeable clay f or maticns. A pressure 
head is created by the slope gradient from the low foothills 
(adyr) along the borders cf the depcesssigns. When these 
horizcns contact sandy gravel deposits, numerous springs 
occur, 

Groundwater of the Ncrthern Eergana depressions 
(Kassansay, Gavasay and ether river basins) is found mainly 
in alluvial and proluvial gravel deposits. In the north, 
cio^e to the foothills cf the adyr, the gravel beds are 
exposed; in the south they ate overlaid by sandy loams and 
loams, Greating donfined aguifer cenditiens. The thickness 
of the Icam depcsits varies from 2-3 m in the north to 20 m 
or more at the Syr Darya. The alluvial fans of the rivers 
in the study area are inccapletely manifested, since the 
rivers have swept away their peripheral parts. Within each 
alluvial fan there is a specific zone of ground and sub-head 
water flow. 

The alluvial fan of the Kassansay Eiver, south of the 
Kassansay and Namanganskiy adyr, is 11-12 km wide and 20-25 ; 
km Icng. Its scutfaern liait is found between the first and 
secend terraces above the flccdplain cf the Syr Darya. Here 
the flow cf groundwater is supplied mainly by infiltration 
from the Kassansayskaya trans-adyr and the Kukuabayskaya 
mezh-adyr depressions. It is additionally reGharged by 
iafiltraticn from surface streams, which form a complex 
hydrcqtaphic retwerk on the alluvial fan. In the upper part 
of the fan, :the water table lies at a depth of 40-60 m; this 
depth decreases as the flew approaches the Syr Darya. 

Water also flows along faults or thoujh the karst 
regiens of carbenate rocks, water . along faults (heedering 
the valley) in addition to the underfloa discharge from 
mountain slopes and valleys raises the level of water bodies 
in the tntee categories discussed above and leads to the 
salirizing and swamping cf iGwlands (Suslov, 1961) . Jcinted 
metam,cr phic a nd igneous recks ha ve stcucturaliy controlled 



Watec pecmeability. They are overlain by shales, with 
interbeds of sanustcne. Grcunduater in these rocks is less 
abuncant tban in the liaestcnes and is characterized by low 
mineral content, monolithic aquifers, and less fissuring to 
provide storage. 

Mineralization of the groundwater confined in 
Mesozoic and Cenozoic fcrnaticns gradually increases with 
distance from the point or recharge. However, the mineral 
content of water emerging at springs does not exceed .5 
q/liter; the water type is calcium hydrocarbonate and 
calciaffl-fflagnesiua hydrocarbonate. 



2 - Adyr area of alluvial terraces; 3 - Moving and fixed sand; 4 - Solonchak; 5 - Marsh; 6 - 
Cultivated oases; 7 - Fergana Canal (from Suslov, 1961). 
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Sncw Hjurclogj 

Btecipitaticc rails as snow loaial^ from Novenbsr 
through March. Cn the souchera slopes of the Coathal 
Raaqe, xroa 25 tc 45 days %ith sncw are observed at the 
1000-1500 m elevations. 

Snow depth and the water content of the sncwpack vary. 
In the middle elevation 2 cne, snow deptn is typically 80-90 
cm, though in some years it ceactes 100-110 cm. On the 
southeastern slopes of the Chatkal Range, at elevations from 
1400 tc 2200 m, the water content of the pads varies from 
100 to 225 mm. In some years in the Padshaata River basin 
it reaches more than 375 om. Figrre 15 shows the 
distribution of water content of snow cover by elevation in 
the GaVasay River in 1942. 

Dates of ablation vary with elevation. In the 
foothills, at elevations tc 1500 m, the snow cover typically 
forms during the seccnd and third ten-day periods in 
November and disappears during the second and third ten-day 
periods of March. At elevations of 1500-2500 a, snow cover 
forms at the end of Cctcber and persists till the end of 
April. At elevations above 3000 m; snow covers the ground 
from the second ten-day period of September until the first 
ten-day period of June. 

Streams in the Syr Darya basin ate primarily supplied 
by meltwater. The indicator of the rivet recharging type 
was calculated as the ratio cf the discharge in July to 
September to the discharge in Match to June. Using this 
methcd, the rivers in this region can be classified as 
rivers whose annual discharge is^ determined mainly by 
fluctuations in precipitation. The effective precipitation 
in the area is retained in storage as sncw and released as 
snow melt periodically to the b>isin's channels (the major 
ones being the l^yr Darya, Naryn, and Sara Darya Riverb). 

The Syr Darya and Naryn have mixed giaciai-snow feeding 
resulting in two high water periods: in spring (April-May) 

from melting snow, and in summer (June-JulyJ during the 
period or ice thaw (Suslcv, 1961). Their tributaries are 
supplied primarily by snowmelt in the spring. In Ap^il 
these rivers manifest violent activizy, yet by July an 
abrupt reduction cf discharge is cbserved due to depletion 
of seasonal snow cover. The variability of annual discharge 
in srow-fed streams is determined primarily by the variation 
in the water content of snow at fhe begiuning of sncwmelt. 
The variability In meltwater froit perpetual snows and 
glaciers, on the ether hand, is net dependent on the ameunt 
of annual pteGipitaticn , but is deteemined largely by the 
amount of energy received Ly the surface of the snow ana 
ice. (Tables 4 and 5) ; 

Because meltwater is the main source or river recharge, 
the distribution of discharge through the year is determined 
mainly by the accumulaticn of sncw and ice in the racantains. 


and partly on variability in infiltration and evaporation. 
Not only the aaplitude of flnctuaticns in the annual 
discharge^ but also the Icny-terji trend of discharges froa 
Central Asia streams are determined by their supply source • 
aeltbater trcn seasccal snou and/cr from perpetual snows and 
glaciers. All the streams iu the Fergana Valley are 
classified as sncv'’glacier and sncK>‘fed streams. In the 
Bountaih re9l<3QS amount ct precipitation generally 
increases Vitb elevation, while vapor pressure decreases 
with increased elevation and lower teaperatures (figure 10). 
For this reason the discharge norm typically increases with 
elevation. However, since the avcunt of precipitation and 
discbarge in the aountains is also influenced by the 
orientation of the slopes in relaticn to the prevailing 
direction of the cove sent of axe masses, the dependence of 
discharge on elevaticn shews up mest clearly for relatively 
restricted otohydtographic re^ichs. 

Table 6 shows the distributuicn of saow surveys in the 
study area. 


o 


o 


o 


Q 



min 



Figure IS. (/Variation of water content of snow cover wi th al ti tude . 
(Basin of Gavasay River, 1942) 

1 = water content at end of February 

2 - water content at end of March 
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RIVER- POST 


Gavasay 
River, 
Mouth of 
the Ters 
River 

Padshaata- 
Mouth of 
the Tostu 
River, 2,6 
km above 
Mouth 


ADOPTED 
GALCULATED 
PERIOD OF - 
OBSERVATIONS 
OF DISCHARGE, 
YEARS 


1937-60 


^tx ^t 


PERIOD 


37 0.33 0.80 


1936-60 


.48 -.45 0.68 0.21 


r ^ 

METEOROLOGICAL 
STATIONS FROM 
WHICH 

METEOROLOGICAL 
FACTORS ARE 
RECORDED 

ELEVATION OF 
METEOROLOGICAL 
STATION, m 

mouth of the 

1818 

Ters Ri ver 


mouth of the 

Tostu River 

1529 


Table 4. Coefficient of Correlation between Discharge of Water and Meteorological Factors 
(Precipitation and Above- Zero Air Temperatures) 
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RIVER 

STATION 

AREA OF 

WAIERSHED 

km"^ 

MtAJI 

ELEVATlOri 

OF 

MATERSUEO 

OBStRVATIOH 

PERIOD 

TOTAL 
NUMBER 
OF YEARS 
OBSERVED 
USED IN 
STUDY 

MEAN DURING OBSERVATION 
PERIOD 

Velocity Modulus layer of 

of Water of Discharge 

flow Discharge nin 

m /sec 1 iters/sec 

per km 

MEAN DURING 

waimmi period 

Velocity lictduTus of layer of 
^of Hater Discharge Discharge 
‘Flow HterS/sec vran 

mvsec per km^ 

CAICUIATED 
VALUE OF 

ADOPTED 

VALUE 

V 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

Padshaata 

Mouth of 
Tosto River 

366 

2830 

1925-62 

37 

6.01 

16.4 

517 

6.03 

16.7 

526 

0.20 

0.21 

Chanacfi 

Chanachsay 

133 

2630 

1938, 

1941, 

1944, 

1949, 

1950, 
1952, 

1955, 

1956, 
1959-62 

12 

(1,76) 

(13.2) 

(416) 

(5.76) 

(13.2) 

(416) 



Ka$sansay 

Mouth of 

1240 

2480 

1945-62 

18 

8.93 

7.20 

227 

(8.50) 

(6.85) 

(216) 

0.26 



the Uryukty 

1780 

2330 

1925-42 

16 

9.80 

5.50 

373 

(10.3) 

(6.07) 

(191) 

0.27 

0.30 


River, 

Baymak 

kishlak 

1780 

2330 

1925-60 

1962 

35 

10.8 

6.07 

191 

lO.B 

6.07 

191 

0.28 

-- 

Uryukty 

Oryuktysay 

86.5 

2180 

1941, 

1947, 

1948, 
1950, 
I960; 
1962 

7 

(0.64) 

(7.40) 

(233) 

(0.62) 

(7.17) 

(226) 



Uryukty 

Mouth of 

Oryuktysay 

River 

116 

2190 

1947-62 


0.48 

4.14 

130 

(0.45) 

(3.88) 

(122) 

0,37 


Alabuka 

Alabukasay 

208 

: 2520 

1938, 

1 1941, 
1944, 
1952, 
1955, 
1957-62 

18 

(2.53) 

(12.2) 

(383) 

(2.40) 

(11.5) 

(363) 



Sunsar 

5Ufnsars.ay 

90>0 

2600 

1941, 

1945-50, 

1959 

8 

(1.01) 

(11.2) 

(353) 

U 

(1.00) 

(11.1) 

(350) 


■ ■ ■ 

Kcksarek 

Koksareksay 

96>2 

2510 

1941 

1 

(0.78) 

(8,11) 

imy 

(0.90) 

(9.36) 

(295) 



&avasay 

Mouth of 
the 

Beshtashata 

River 

288 

2840 

1932-38 

6 

(3.82) 

13.3 

419 

(4.20) 

(14.6) 

(460) 



Gevasay 

Mouth of 
the Ters 
River 

361 

2840 

1932-62 

31 

4.96 

13.7 

432 

5.05 

14.0 

441 

0.28 

0.28 

Gc vasay 

Gava 

kishlak 

657 

2460 

1925-62 

36 

5-93 

9 10 

287 

6.29 

9.57 

301 

0.31 

0.32 
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General Hydrclo«jy 

Streams in the study area ate also supplied, tc seme 
extent by inflow ox gtoundHater from mountain valleys and 
slopes. This scuEce is also indirectly dependent cn 
snowmelt, since groundwater recharge is uetermineu priaatxly 
by the vclume of meltwater. Streams fed from groundwater 
occur in the upper foothill belt (up to 1000 m) and are fed 
by groundwater and showers. They tend to have small steady 
or intermittent discharges with cccasional flash flooding 
during intense storms (Suslov, 1S61). 

1 The Syr Carya follows the nerthern rim of the Fergana 
Villey. Its upper section, the Karyn Biver, enters the 
valley from the east through the Fergana Range. Surface 
water in the basin cccuts as flow in the Syr Darya, in its 
tributaries (including the Sexh River), in the dense netwcrlc 
of itrigaticn, and as standing water in shallow closed 
basins and artificial reservoirs (Figure 2) . 

In the study area, the Gavasay, Surasar, Kassansay and 
Padshaata rivers have the highest discharge of streams which 
drain the southeastern slopes of the Chatkal range and the 
Kuramin mountains. All these rivers are right banlc 
tributaries of the Syr Darya. These rivers do not freeze 
over widely, (the period when these streams ace frozen 
typically does net last longer than 30-90 days) , hewevee, 
local freezing car lead to the fermation of major ice jams 
which create significant fluctuations in water level daring 
their formation and breakup. Foriaticn cf ficating ice is 
prometed by high flew velccities, cocky, ragids-filled 
water, etc. These ice focmatiens can last two or even three 
months . 

The turbidity of water in the Chatjcal Range varies from 
less than 0.1 kg/ni3 to 0.3 kg/m^. In spring and the first 
half cf summer, all rivers draining these mountains carry 
mudflows, often with destructive force. 

Figure 16 shews that, at the same elevations, the 
Padshaata River basin is characterized by greater discharge 
moduli than the basins of the Gavasay and Kassansay Rivers. 

At the elevation cf 1600 m the discharge modulus of the 

River is 1.0 litec/sec while in the Gavasay 

and Kassansay rivers it is 0.5 liters/sec km^. This is due 
to the areal distributicn cf precipitation, which varies 
with the orientation of the basins. The discharge cf the 
Gavasay and Kassansay rivers, whese valleys are Unfavorably 
oriented, averages S liters/sec km^. Further to the cast, 
ociertaticn is mere favorable , so the moduli of the 
discharge of the Padshaata River and several left lank 
tributaries of the Kassansay Rivet increase to about 1 1 
liters/sec koi^. 

The period of greatest discharge frcic streams in the 
study area cccurs frctt 2aich to August oc September. 75-35^ 
of the anaual discharge occucs at this time. The highest 
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aonual flow velocities typically occur in ilay; however, the 
discharge of the Eadshaata fiiver is highest in June. Low 
water lasts fron Septeaber tc March. Coefficients cf 
variation, describing the annual distributuion of flow, are 
shown in Tablh 7. Table 8 shows the aonthly distributuion 
of discharge as a percent cf annual discharge. The 
relatively low elevaticns cf the Chatkal and Kurattio Ganges 
account for the similarity in the discharge of the stceaas 
in the study area. Generally, then, three periods can be 
distinguished: 

1. The period of sncw high water, formed mostly by the 
meltwater from seascnal sncws at low to mid-elevations. The 
onset of this high-water period is associated with the cnset 
of stable, above-sere tempecatures. The volume of discharge 
is determined by the amount of snowfall of the given year. 
This pecied typically lasts from Harch-June. 

2. The period of sncw-glacier high water, formed 
mainly by^ the meltwater from high-acuntain snows and 
qlaciecs. This period coincides with the hottest part of 
the year (July- Septs abac) . 

j! 3. The period cf low-water, when the streams are 
supplied primarily by greundwatet. This period is 
characterized by stable, relatively low discharge and by the 
absence of diurnal flactuaticns. It typically lasts from 
Oct cber-Fetruary . 

The characteristics cf high-water and low-water periods 
Of the rivets in this- region, and the maximum and ainimum 
volume flows ate in Tables 9, 10 and 11. Table 12 shows the 
types of measurements taken at the hy dcoieteorological 
stations in the study area. 

Movement of surface watec and groundwater in the 
Fergana Basin may be described by an inf low-outflcw budget 
acccunting for the differert modes of water transport. The 
primary inflow compcDents|precipitation, streamflcw and 
qroujidwater) have been described above. The total water 
resources of the entire Syr Carya basin in 196d were 
estimated to be about 35 sq km (Keitskiy, 1968) . About 
two-thirds cf tbe arriving runoff is expended on irrigation 
of fields, evaporation from the river Ghannels and 
reservoirs , and the periodic flooding of lowlands. Mater 
loss also cccurs by infiltration (undecflcw or thrcughflQW) 
evaperatien from surface water and evapotraDspiraticn. The 
remaining discharge flews to the firal sea.; 

The main inflow compenents of the Fergana Basin may 
then be summarized as: I) river flew £rani tne;1{ara Darya 

and Natyn, 2) stream and river flew fron the surrounding 
mountains (Tien Shan, Eamir- Alai) , 3) effective 
precipitat ion in the ierm cf snow and ice melt, and 4) 
undetflow cr thrcughflcw cf greundwater froai tae m<|untains. 
The main outflow components include; 1) discnacge, cf the Syr 



Dacya surface waters out cf the Fergana Basin, 2) 
infiltration or underflow losses cf groundwater out of the 
basinv 3) evaporaticr fros standing and flowing surface 
water, and 4) agricultural and vegetation lasses 
{evapottanspiraticn) » ^ 

Careful nanageaient of water resources in the Fergana 
Basin is required to oarinize the rescurce benefits. An 
increase in ccnsumpticn of water can promote conditions in 
which irrigation needs arei'nct met. This is shown by the 
fact'^tbat the water level in the rivers does not always 
eniute the entry of the water at the heads cf the irrigation 
systems. A considerable pcrticn of the runoff in the area 
passes in the form of high, brief high-water peaks 
superimposed on the meltwater wave. Cn rivers which are fed 
by runoff from relatively low mountains or small drainage 
basins, much water passes at the beginning of the warm 
season , prior to the basic irrigation period. 


1.0 


1.8 


2.5 


3 4 H 


CP km 


Figure 16. Dependence of modulus of discharge on mean elevation of 
watershed (km) 


A. Padshaata River 

B. Kassansay River 

C. Gavasay River 



RIVER. STATION 

Padshaata, mouth of the Tostu River 
Kassansay , Baymak kishlak 
Gavasay, mouth of the Ters River 
Gavasay, Gava Kishlak 


COEFFIC 

Spring 

0.23 

0.32 

0.28 

0.35 


OF VARIATIOM OF SEASONAL DISCHARGE 


Summer 

Autumn-Winter 

0.27 

0.22 

0.52 

0.29 

0.44 

0,29 

0,44 

0,26 




RIVER, 
STATION, 
WATER LEVEL 


OF YEAR 

Jan 

Feb 

Mar 

Apr 

Padshaata, 
Mouth of the 
Tostu River 
High Water 

2 A 

2.1 

o ^ 

2.7 

8.7 

Medium 

3.3 

2.9 

3.1 

7.8 

Low Water 

3.1 

2.8 

3.4 

; 6.8 

Very Low 
Water 

3.3 

2.9 

3.4 

6.8 

Kassansay 
Bayinak 
kishlak 
High Water 

2.7 

2.2 

1.9 

6.7 

Medium 

2.8 

2.7 

3.1 

6.1 

Low Water 

3.0 

2.7 

4.0 

9.0 

Very Low 
Water 

3.1 

2.8 

4.6 

10.2 

Gavasay, 
Mouth of the 
Ters River 
High Water 

1.4 

1.4 

1.8 

9.8 

Medium 

1.9 

2.2 

2.2 

11.6 

Low Water 

1.9 

1.7 

3.3 

9.6 

Very Low 
Water 

2.1 

1.9 

3.4 

9.7 


M O N T II L Y D I S C \\ A R 6 E 


May 

June 

July 

Aiig 

Sep 

Oct 

Nov 

Dec 

17.0 

21.9 

17.8 

10.5 

5.8 

4.1 

3.8 

3.2 

17.8 

21.8 

16,8 

10.2 

6.0 

4.0 

3.6 

2.8 

17.1 

23.5 

15,8 

9.6 

6.0 

4.6 

3.7 

3.6 

17.3 

23.7 

15.1 

9.2 

5.7 

4,8 

4.0 

3.8 


15.9 

22.7 

18.1 

14.9 

6.4 

2.9 

2.4 

3. 

17,4 

23.2 

17.0 • 

11.9 

5.2 

3,7 

3.7 

3. 

16,9 

22.1 

14.5 

8.6 

5.6 

5.J" 

5.0 

3. 

19.2 

25.1 

10.6 

6.3 


'^l4 

5,2 

3, 









33.3 

25,8 

11.9 

4.9 

2.8 

2.6 

2.4 

1, 

23.3 

34.6 

9.9 

4,6 

3.0 

2,5 

2.5 

1. 

35.9 

23,5 

8.2 

4.3 

3.1 

3,2 

2.8 

2. 

36.2 

23.7 

7.2 

3.8 

2,7 

3.4 

3.2 

2. 
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RIVER, 
STATION, 
WATER LEVEL 
OF YEAR 

Jan 

Feb 

Gavasay, 

Gava kishlak 
High Water 

1.7 

1.6 

Medium 

2.1 

2.0 

Low Water 

2.5 

2.3 

Very Low 
Water 

2.9 

2.6 


Chaadak 
(Chaadaksay), 
Mouth of the 
Dzhul aysay 
River, 


High Water 

1.5 

1.7 

Medium 

1.9 

2.0 

Low Water 

2.3 

2.7 

Very Low 
Water 

2.9 

3.5 


MONTHLY 


Mar 

Apr 

May 

June 

2.2 

12.2 

30.9 

23.9 

2.5 

10.3 

32.4 

23,7 

3,6 

8.8 

32,8 

23.5 

3.6 

8.8 

32.6 

23.4 


3.2 

15.6 

24.6 

31.3 

3.8 

14,0 

33.6 

22.6 

4.8 

11,9 

36.6 

20,3 

4.7 

11.7 

35.8^ 

19,9 


O 


0 




DISCHARGE 

July Aug Sep Oct Nov Dec 


11.5 

5,3 

3,0 

2.7 

2.8 

2.2 

9,9 

5,1 

3,5 

3.0 

2,9 

2.6 

8.9 

4,4 

3.4 

3,5 

3.4 

2,9 

8,0 

4,0 

3,0 

3.9 

3.9 

3,3 


9,35i 

4,1 

2.3 

2.2 

2.3 

1, 

8.0 

3,9 

2.7 

2,8 

2.5 

2. 

6,8 

3,4 

2,8 

2.9 

3.0 

2. 


5.6 2,8 


2.3 3.7 3.9 3,2 


CT> CM in 









RIVER, 

STATION 


Mouth of 
the Tost! 
River 


Mouth of 
the 

Uryukty 

River 


Baymak 

kishlak 


Gavasay, 
Mouth of 
the Ters 
River 


PERIOD OF 
OBSERVATIONS 

CHARACTERISTICS 

"Years Number 


of 


Years 

u;- ■ . 

, 1926-62 37 

^max 


^mt 


h 

, 1952-62 11 

^max 


*^mt 


h 

, 1926-41 16 

^max 


^mt 


h 

1933-62 30 

^max 


^mt 


h 


DURIfIG PERIOD OF OBSERVATION 
Highest Mean-Hultiannual 


3 

Q, nr/ sec; 


3 

Q, ni /sec; 

Modulus of 

h, mm 

Year 

h» mm 

Discharge, 
liters/ sec 

44.5 

1928 

28,1 

26.8 



22.5 

61.4 

623 

1960 

408 

9m 

84.8 

1953 

51.6 

41.6 

-- 

-- 

46.7 

37.6 

161 

1952 

104 

— 

132 

1928 

68.3 

38.4 

239 

1934 

132’ 


(55.7) 

1934 

38.0 

105 

— 

— 

29.8 

82.4 

601 

1949 

362 

— 



0.2a 

0.27 

0.24 

0.30 

0.30 

0.^36 

0.40 

0.36 


&S 

as 


Table 9. Maximum Annual 
FI o\i Rates of,Meltwater 


Velocities of Water Flow/ (Q . , m /sec). Maximum Discharge 
(Q 4 ., m^/sec) and Discharge Piper (h, mm) in Freshet Period 
















RIVER STATIOM AREA OF PERIOD OF 

WAIERSHEO OBSERVATIONS 


NUMBER OF LOU-HATER PERIOD 

OBSERVATION Mean Date Duration* Mean Layer of 

YEARS Beginning Ending Days Consumption Discharge, 

* During Low mm 

wf 

\^rP«r1od 

■ . '' mVses={. 


Gavasay 


Gavasay 


Mouth 

366 

1931-34 

of the 


1936-48 

Totsu 
River, 
2.6 km 
above 
mouth 

I' . ■ ; ' 

1951-62 

Mouth 

361 

1936.19 • 

of the 


1939-62 

Ters 

River 



Mouth 
of the 
Gava 
River 

657 

1936-62 


3/24 190 


3/24 221 


3/27 209 


Table 10. CharactefistiGS of Low Water Period 




Gav^say 


Gavasay 


ve 

the mouth 


Gave 

Mshlak 


MEAN 

ELEVATION 

OF 

WATERSHED 

m 

PERIOD OF 

OBSERVATIONS 

(ASSUMED) 

NUMBER 

OF 

YEARS 

MINIHUH 
DAILY 
DISCHARGE 
1 Iters/sec 

2830 

1936-62 

30 

4.46 

2840 

3937-65 

29 

1.88 

2460 

1936-66 

31 

1.45 


MONTHLY MEANS FOR PERIOD 
OF OBSERVATION 
Consumption Modulus of ^ 


ot water 
m /sec 


Discharge 

Lijers/sec 


COEFFICIENT 

OF 

VARIATION 


Table 11. Mean Mon Diurnal Discharge 
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RIVER STATION 

OISTAliCE: 

AREA OF 

PERIOD OF OBSERVATION 

MAXIMUM 
vri nriT i rc 

PRESENCE OF 
Mr Tr noni nr ttai 



fRON 

MOUTH, KM 

UAJERSHED 

Level 

Velocity 

Temperature 

Glaciation 

Suspended 

Detritus 

Chemical 

Composition 

TCt-UL 1 1 11.0 

OF HATER 

now 6 

LAYER OF 
fRESMET 

nt 1 tUr.ULUUi LnL 

STATIONS (M/CT) 
AND POSTS (H/P) 


_J 














3 

4 

5 

6 

7 


9 

-IQ „ . 

11 ; 


Radsliaata, nouth 


^66 

1^54-64 


1953-76 

1936-69 



1926-33 

M/CT 


of the Tostil 









1935-48 


River, ?/6 km 









1950-70 


■■ 

above moutit 












Padshaata, 

AZ 



_« 

1941-51 

■ 



-- 




Yangikurgan 

kishlalc 




1953-62 








Chanach, 

15 

133 


1938-41 


-- 

-- 

-- 

-- 

-- 


Chanachsay. 8 km 




44,49.50 








from the village 




52,55,56 








af Aktam 




59.68,70 








Kassansay, 
Ky^yi- Tokay 
kishlak 

74 

1130 

, ■ 

1967-70 

— 

1967-70 

1968-70 

1968-70 

1967-70 

— 


Kassansay, 

66 

1240 

1951-67 

1945-68 

1952-68 

1951-68 

1956-68 

1962-68 
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Table 12. Extent of Hydrometeorological Coverage 
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Irrigaticn and Agriculture 

Approximately 27% of the total laad area of the USSR is 
agricultural/ and slightly mere than cne-tiiird cf the 
agricultural land is araisle. The remainder of the 
agricultural land is meadc!i</ pasture/ orchard/ vineyard, or 
is idle. In 1972, cf the 224 million hectares of arable 
land in the USSR, almost all Has seen. Irrigated land at 
this time totaled 12 millicn hectares, nith 43% 

(approximately 5 million hectares) cf this located in 
Central Asia, Hhete irrigaticn has played a major rcle in 
agriculture fer centuries. As of 1970, over 800,000 
hectares were irrigated in the Fergana Valley, representing 
more than 95% cf the arable land (Lydolph, 1970). In the 
valley much of the water from the mcuntain streams is 
absorbed by the loose material of the alluvial fans, and the 
sediments deposited result in the formation of fertile soils 
of great agricultural value- Intensive cultivaticn by a 
large population has always been possible on the oases found 
at the foot of the large alluvial fans in the valley. The 
natural irrigation system of the alluvial fans has recently 
been improved with the ccnstructicn of numerous irrigation 
canals. Because cf the increase in irrigated land in the 
last 35 years, planting is new dene on mountain slopes, 
along the strips bordering the canal paths, and in the 
glades between orchards. 

The northern part cf the Fergana Valley has favorable 
n^^atural conditiers fer intensive cultivation and irrigated 
agriculture. The main methods of irrigation are furrow fer 
plowed crops, and flooding along strips for cereal crops and 
alfalfa. The irrigated fields ate tilled entirely 
mechanically. 

Agriculture in Soviet Central Asia is largely dependent 
on : irrigation because evaporation exceeds precipitaticn over 
most of the area. Snow and rain in the mountains to the 
south feed the large rivers that flow into the lowlands, 
depositing alluvium and providing water for extensive canal 
netwerks (CIA , 1974) . 

The main canal system ccmpletely encircles the basin 
(Figure 17) , with feeder canals distributing water over much i 
of the valley fleer. Only the lowest central pertien of the 
basin is not irrigated due to texture and brine conditions 
unsuitable for agriculture. Much cf the water supplied to 
this part of the basin occurs in the form of groundwater. 

Almost all flew from the fcothill tributaries on the margins 
of the Fergana Valley is diverted into canals for 
distributicn as irrigation. , Little, if any, reaches the Syr 
Darya. 

The largest canal is the 30C km Great Fergana Canal, 
fed by the Karyn River and Kara Cara. It runs the length cf 
the valley from the northeast to the, southwest and fortes the 
base lor the irrigation system in the Fergana Basin (Figure 
17). The Forgana Canal passes through the zone of laximuir 
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qcouDdwate'r, skirting the ferioetec of the allavial fans in 
the valley. & oumbei cf hydrcelectrical iastallaticns, 
locks, sluices and other uater works Ci:erate along its 
length (Wheeler and Fcctnan, 1954). The function of the 
Canal and its subchannels is to sufplement pre-existing 
irrication systeas on the northern slopes of the Alai Range 
tith water cf the Haryn River and Kara Darya. The sain 
canal intersects all of the irrigated fans ox the rivers and 
streams in the southern part of the valley. The canal 
system begins with water from the Haryn which is transferred 
along the canal tc the Kara Darya, dammed by Kugan-’iarskaya 
Dam, which raises its level by 4 m to allow water acvement 
under gravity westward. ^^Dnring the period 1957 to 1959, a 
large irrigation dam was\ccnstructed near the apex cf the 
Kokand alluvial fan to collect a portion of the waters of 
the £okh River for the ifrigaticn of an additional 60,000 
hectares or land in the ^cuth-central section of the valley. 
In addition tc water from the Haryn and Kara Darya, the 
canal system receives additional water from rivers flowing 
off the mountains of the Alai Ranee, and minor amounts from 
flow from the Tien Shan. Rivers cf different regimes are 
united by the system. Rivers cf predominant glacial feeding 
(Sokh, Isfara) with the main high water in summer, rivers of 
snow feeding (Kara Darya) with main high water periods in 
spiring, and rivets of mixed feeding (Haryn, Syr Darya) with 
bimodal high water periods are joined in a surface water 
system with an almost completely adjusted regime with 
relatively cerstant flew spring and summer (Suslov, 1961) . 
Many feeder canals distribute the water ever much of the 
basin floor (S5X cf arable land) (Gregory, 1968), 

The Northern and Southetn Fergana Canals are the two 
other major irrigation systems in the valley. The South 
Canal is about 105 ko long and irrigates the land along the 
foot of the meuntains in the southeastern pertien of the 
valley. In addition, many streams hhich flow into the 
valley and have built alluvial fans are used for gravity 
flow irriqaticn. 

The area included in the Soviet study area is served by 
several irrigation systems; the Northern Fergana Canal, 
which irrigates 63,000 hectares; the Radshaatinskaya Canal, 
which irrigates 32,000 hectares, the Kassansayskaya Canal, 
whicb irrigates 24,000 hectares has a reservoir with a 
165 million m^ capacity; and the Gavasayskaya canal, 
serving 17,00C hectares. 

The North Canal, 162 km long, pcovides water to the 
area north of and parallel tc the Syr Darya. It has existed 
in ixs present form since 1960, and received water from the 
Ochkurganskaya storage dam on the Hafy^n siuce 1967, 

The irrigatioE system serves a strip from 2-6 km wide along 
the canal. Tbe canal is built along a hillside and 
traverses numerous ravines and streams. In winter the water 
supply is hampered by the fetmatien of ice») ^ 
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Ihe Padshaatinska ja sjstea is fed bj the Padshaata 
Bivec. Springs are an imjfcrtant source of water to this 
system* The total discharge ox the springs is about 160 
million m^; the mean discharge of the riwer at its 
emergence from the mcuntains is 190 million There are 

also 100 wells with a tctal jield of 2-3 m^/sec- The total 
length of the canals in the network is about 1100 km, of 
which 160 km make up the large trunk and interfars canals. 
120 km of these are concrete- lined, 

■ The Kassanskaya irrigation system is fed by the 
Kassansay River « whose source is above 4000 m. Its 
reservoir allows seasonal regulation of flow. The trunk and 
interfarm canals are 176 km long, of which about 40 percent 
are concrete'- lined. 

The Gavasayskaya system gets water from the Gavasay 
River « which is fed by seascnal snows. Springs and wells 
are important to the system - springs supply about 1.5 
mVsecj around 100 wells have a total yield of 3 m^/sec. 

Of the 97 km of interfarm canals, 45 km have concrete 
linicg. 

Cultivation and processing of cotton, the valley’s most 
valuable crop, are the chief economic activities in the 
Fergana Valley, and Kokand, Hargelan, and Fergana are 
important cotton-ginning centers {Plummer, et al., 1971). 
This basin produces between two-thirds and three-guarters of 
the cotton grown in the Soviet Onicn {Gregory, 1968). 

Alfalfa is grown for livestock fodder and to reduce the 
inherent salinity of the soil. Wheat, barley, corn, 
sorghum, melons, and vegetables ate also cultivated. 
Vegetables are grown fcr local urban markets. Other food 
grains are cultivated without irrigation. 

The Fergana Valley is one of the main orchard- vineyard 
districts cf Central Asia. Fruit crops, such as apricots, 
peaches, and plums, are prcdoced, particularly around the 
larger cities. Fruit is also dried in local factories 
{Gregory, 1968) . The climate of the Fergana is very 
favorable for the cultivation cf grapes {BuriJtin, 1S67). 

Rice is another major crop of the valley, although a minor 
food crop for the Soviet Union as a whole. Sxlk is produced 
at mulberry plantations in the valley (Figure 18). 

Livestock prcducticn in the Fergana Valley includes 
both beef and dairy cattle breeding, and sheep breeding for 
muttcE and wool. Ihe distribution cf meat-producing ani.ma3.s 
in the USSR follows the distrihuticn of population, unlike 
the situation in the United States, where livestock tend to 
De concentrated in less populous states;. Dairy cows in the 
USSR are utilized fcr both beef and milk, and the dairy 
industry is population-oriented. ihe Soviet livestock 
industry has alsc been labor-intensive and thus is well 
suited to areas cf dense rural settlement {CIA, 1974), such 
as the Fergana Valley. 
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Industry and Manufacturing 


The Fergana Valley ccnstitutes an iapoctaut econoaic 
district of the Uzbek, TadzhiJc and Kirgiz SSfi's. The 
coobinaticQ o£ coal and iron deposits has contributed to the 
area's industrial developaent. dicng the industries in the 
area are oil, aining, engineering, chenicais, textiles 
(cotton and silJc) ,. food processing (including vegetable oil 
and vine) , and fertilizers. In spite of the Fergana's 
industrial productivity, Uzbekistan is not self-sufficient 
in steel, fertilizer, or fuel production (Gregory, 1968). 
Much of the industry and oanufacturing in the Fergana Valley 
is located in a ring around the valley floor which follows 
the Bain roads and railroads that skirt the alluvial fans. 
Industry tends to be in greater concentration along the 
soathern portion, and especially at the eastern end in the 
Andizhan region. The Kokand, Nasangau, and ilargelao-Fergana 
regicns also include industrial ccocentraticns. 

Industrially, Andizhan is the aost diverse city in this 
region, being involved with the: food; building materials; 

53 heoical; agricultural engineering; engineering; 
and setal- working industries. Next is Fergana, with food 
industry, silk, textiles, chemical industry, and oil 
refirging. Kokand has food industry and textiles, chemical 
industry and engineering and metal working. Kokand lost its 
piace as an acministiative center, after the revolution, to 
the rising industrial center of Fergana, and conseguently 
sulffered a drop in population. However, the growth of 
Kokand's cotton-ginning industry and the addition of aetal 
working and chemical fertilizer industries (phosphate and 
super-phosphate) have helped to revitalize the city. 

Namangan is involved with the food industry, textiles, and 
building materials (Burikin, 1967). The food industry, 
although found in most of the major cities, is principally 
concentrated at the eastern end of the valley. 

The cities of Pap, in the northwest, and Shor-Su, in 
th^ southwest, are both involved with the chemical industry. 
KuVasai, in tbe south, has a cement works and building 
materials industry {Burikin, 1967 ) . Begovat was established 
by tie Soviets as the only steel production center in 
Central Asia. In 1968 the annual production of steel was 
about 334,000 tens (Gregory , 1 968) • It is hoped that 
eventually enough steel can be produced to make the Central 
Asian fegicn self— sufficient. 

Fergana Valley centains extensive deposits of oil 
and natural gas, principally along its southern rim. Here 
there are several oil and natural gas extraction centerlE. ! 
Several natural gas pipelines run aiong the southern edge: of 
the valley (Buriki n , 1367) , In tne Sura man Han ge to the 
northwest of the valley fleer are large depoists cf lew 
qrade coal, as well as some copper, dluminam, and other 
minetais in lesser amounts. ’While the ainerai depcsits of 

the Fergana Valley ate net large cn a global scale, they 
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ibstract 


i€ develop a digital conputer nodel, based upon 
eqnilibciun tenpecature theory, o£ the' aagnitudes of the 

I 

conpcnents of the energy balance for a aountainous 
topographic surface. The nodel includes an algorithm for 
calculating slope, exposure, and horizon information for 
eivery point on a digitized grid. Incoming solar and thermal, 
dadiation are then simulated for every point. At present the 
model is designed for clear-sky conditions, but it could be 
modified for cloudy skies. Air tenpecature and humidity 

I 

Variations are specified b^ externally defined relations, and 
values for albedo and soil thermal properties are specified 
for every point on the grid. Wind speed variation over the 
grid is not modeled, but is specified by an empirical 
function, the model simulates net radiation, soil heat flou, 
sensible and latent heat flov, and surface temperature at 
specified time intecfals. Cutput is in the form of contour 
maps. 


Key Hords: energy balance clifflatolcgy, computer simulation. 
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Surface energy balance nodels have been used to sinulate 
accurately the surface temperature, soil heat and water flow, 
ahd snowmelt. Their operational utility over wide areas of 
variable terrain has been limited, however, because of the 
extensive data needed to drive such models. In this paper we 
present a method and a description of a computer prograa for 
extending energy balance simulations over rugged terrain, and 
illustrate a few cases of such siaulations. 
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Fourier sine coefficient 
mass heat capacity 

heat capacity of air at constant pressure 

dust and haze cootent in units of 800 particles cm~*a 

diabatic exchange coefficient 

vapor pressure 

azimuth of slope from south (negative when east of 
south) ' 

Fourier modulus, = a At (As)-* 

gralj^taticsal ‘ acceler atldn 


G soil heat flew 


h gtid spacing 
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B sensible heat flow 

I 

incoming atnospberic (infraccd) cadiation 
j subscript £qc node 

If outgoing thersal radiation 

J Julian date 

k , fon Karnan constant, = 0*4, 

' I ■ , , 

It^. absorption coefficient 

scattering coefficient 
K thermal conductivity 

L latent heat of vaporization, =* 

LE latent heat flux 

a^ molecular Height ^ 

H optical air mass 

Pi^ precipitahle water vapor 

P atmospheric pressure 

g incoming diffuse solar radiation (including 

hackscatter) 

C incoming beam solar radiation 

Q* incoming beam solar radiation on slope 

Cq solar constant » 

r earth's radios vector - f' 

i gas constant 

H„ net radiation, * Q(Wq) ♦ g(1-aq) ♦ 1+ - l+ 

s (subscript) surface 

s depth 

S slope from horizontal 

. J ^ , ... . . . . _ ....... 

t ■ time 

air temperature 
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teapecature. at altitude z 
T^i aeasuced teapecatuce at altitude 
u uind speed 

u« ftictioo velocity 

i ■ 

Vj; tberaal cadiaticn vieu factor 

z altitude of surface or of teapecature, uindr and 
huaidity values 

i 

z^ roughness length 

Z solar zenith angle 

Z* solar zenith angle with respect to slope 
Z|^ horizon angle frca zenith 
a theraal dif fusivity, = K/pC = f (T) 

5 solar declinaticn 

ejj effective atmospheric eaissivity 

Eg surface enissivity 

6 potential temperature 

6^ scaling temperature 

p density 

a Stefan-BoltzBann coDStant 

9 latitude 

0 ) hour angle of the sun from nccn (negative in morning) 

BBVIEH OF SURFACl ENEBGT EALAHCF MODELS AND NUHEBICAL METHODS 

Broadly defined, surface energy oalance models are 
those uhicb simulate all of the major surface energy 
exchange processes, and vhich use the conservation of energy 
principle as a boundary ccnditicn. In its simplest form, 
over a solid surface, the principle is described by the 
equation: 

a„ + a ♦ LE G = 
n 


0 


( 1 ) 


In the soil, heat floii is described by the one* 
dieensional conduction equation: 

3T/at = a£l#s] 3*Vas* 

f ' 



Typically it is assuned that no diurnal heat flux occurs 
belou a certain depth, e.g. 30 ca* 


Ihe sensible and latent heat fluxes are estiaated frqa 
air teaperature and huaidity, surface teaperature and 
huaidity, wind speed, and roughness length, using the 
Businger*-Dyec foraulation for ataospheric teaperature and 
vind speed profiles: 

H * Pa Cp k u* e* {3) 

LE ^ i .622 a L (63-65) (6^-65)-* (4a) 

Xf » 6., the latest heat exchange is 

o 5 

; ^1 ■ 

IE « .622 L X* u («a~®s^ £ln( 2 /ZQ)]-« (4b) 

The friction velocity u* and the scaling temperature 6* are 
givee by 

K . : ' - ■ ■■ '■ ■ ■ . 

t u* » Xu/£ln(z/2Q)-,j,^J V (5) 

e* * (6) 

j 

and \p 2 i Busingec-Cyer \{)-f unctions; their 

formulation and derivatico are described by Paulson 

; 'i ' ' 

only mean values of teapeiatuce, Kind speed, and huaidity ace 
available, the best "calculation height" for z is given by 
equating the bulk adiabatic' diffusivity £29] and the daily 


distocbaace penetration depth diffusivity £46J: 


hgu , (7) 

£ln(z/ 2 Q )32 . 5.184x10« 

This can be solved iteratively for z. Although the height 
estisated by this egoation is typically lower than the 
height of the diarnal heat pulse, it provides an estiaate 
which is consistent with the assuaption that heat fluxes are 
constant with height in the boundary layer. 

In order to evaluate aagnitudes of the various beat 
transfer processes, a powerful solution method based upon 
equilibrium temperature theory has been developed £29]. For 
a given set of astronomical, temporal, atmospheric, and 
surface boundary conditicns, there is a unigue surface 
temperature which will satisfy eguations (1) through (7). 

Halstead et al. £18] designed an analog computer to 
estimate the surface temperature and the other components, 
istogue £12] adapted a similar model, driven by measured 
values of incoming solar radiation, air temperature at screen 
height, wind speed, and sub-surface soil temperature, to a 
digital computer. Hyzup £29] simulated the urban heat island 
with an analog computer, illustrating the iapoL’tance of soil 
moisture in deteraining the disposition of the energy balahce 
and therefore the surface temperature. Sasaoori £44] 
simulated atmospheric boundary layer and soil temperature 
profiles, plus wind velocities in the boundary layer, using a 
digital model driven by measured values of solar radiation, 
air temperature, and wind speed, at a given level. Cutcalt 
£32] devised a digital model which simulates solar and x)Mc 
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infrared radiation, and which nay be driven by eean daily 
tenperatures if hlgher**f repuency observations are not 
available; this coaputer prograa has been revised to use the 
Bjusinger-Dyer functions for turbulent heat and water vapor 
transfer, and adapted to a topographic surface for this 
paper* 

lo solve for the egullibriua surface teaperature an 
4nitial value of is chosen, and values of soil, sensible, 

an^ latent beat flow are calculated. These are used in 

/ 

/ ' 

equation (1), and the process is repeated, with iterative 
variations in until eguaticn (1) is satisfied within a 
specified aargin, for exaaple 1 ia'*^ oc 1 aly ain~>. The 
soil temperature profile is then adjusted by eguaticn (2) , 
coapletlng the tiae step. In the current version of the 
coaputer program we use a secant algorithm £ 1 ] to find the 
equilibrium surface temperature; generally it converges 
within three iterations. The structure of the secant 

algoritha is identical to that of Newton's method, except 

■ I ' . ' ■ 

that the analytic expression for the derivative is replaced 

by; a first-order forward difference approximation. If the 
l.ft ,i4€ of e,u.tio. (1) is O.rtossod Jeoeraily as f(T^». 

then subseguent guesses for T are given by: 

s. 

Ml M*^ ^ £ i IF ^ \ » fft ^ fV n * X ^ ^ t m ^*X.« • 

Ihere exist many schemes for numerical solution of the 
heat flow eguation (2) . The present coaputer program eoploys^^^^^^^^^^ 
a fully implicit method. This has the advantage of 
unconditional stability, regardless of the value of the ^ 



Fourier nodulus F ~ a At (AS)"^* and the choices o£ At and as 
are governed only by the accuracy desired £423* The soil 
thernal diffusivity a is a function of teaperature and eater 
content £2* 47] and thus varies with depth and tiae« The 
iaplicit scheae used is £42]: 


where subscript j indicates the node and superscript t 
indicates the tiae step. Application to all nodes at one 
tine step leads to the following tridiagonal systea of 
siaultaneous linear eguations: 


0 0 of' 


fr 



-Fg^ U2F^^ -Pg* 0 0 


T t+1 
■^2 



0 -Fjt uaij* -?3‘ 0 


T t + 1 

^3 



0 0 . • . 


• 


• 

0 0 0 .. 


• 


^ * > 



Teaperature at the top node is deterained by eguation (8) and 
teaperature at the bottom node is held constant. 

Hany experiments have demonstrated that such energy 
balance models simulate surface temperature with reasonable 
accuracy, in urban £29, 33], mid-latitude £44j, arctic £34, 
35], and alpine £5] envirccments. These simulations have 
been verified by point surface temperature measurenents and, 
in the case of £33], by thermal imagery. 


Such models have also been coupled with simulations of 
soil water movement to produce combined heat and water flow 
models [ 17, 36, 43], These have importaut iaplications for 


coDStruction lo axctic areas, uhece the noveoent, freezing, 
and thawing of soil water nay have severe conseguences. 

Energy balance sinulaticn has also been applied to 
developoent of accurate snowaelt aodels. Outcalt et al. [37] 
used station weather data to siaulate changes in the snow and 
soil thermal regimes at Barrow, Alaska. They simulated both 

i 

the beginning and end of the melt period within cne day of 
the actual tiaesr and the depth of thaw penetration into the 
sbil to within 5 cm. Anderson [3] compared the results of 

hiis snowmelt simulation with measurements of snow cover water 

i 

equivalent at the National Oceanic and Atmospheric 
Administration * Agricultural Besearch Service cooperative 
sinow research station near Danville, Vermont, and found that 
in most situations his energy balance model was superior to 
more conventional sncvmelt aodels. 

In cases where surface temperature is known, for example 
from thermal imagery from spacecraft or aircraft, energy 
halance models car be applied in reverse to estimate the 
thermal characteristics of the surface material, which can 
then be used to identify the surface lithology [21, 22, 39, 
46]. Alternatively, if "the surface lithology is also known, 
thermal imagery can be used tc locate areas of high 
gebthersal heat flow £9, 48}. leSchack and Del Grande £26] 

■ ' I , 

i 

have outlined a method whereby thersal sensing in two 

i 

wavelength bands (e.g. 5 and 10 ym) can separate 

! ' 

enlissivity-related effects fees tespecatuce> related effects, 
thereby allowing determination of actual temperatures to an 
accuracy of 0. 1 or 0.2 ®K. 
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Ihe application of an energy balance node! to a 
topogiaphic surface necessarily involves the use of a 
digital conputcr* Therefore the tern "topographic surface"^ 
as used hence, refers to a digitized grid, consisting of a 
finite nunber of points- Such grids are now available at 
very nodest cost as "Digital Terrain Tapes" fron the 
National Cartographic Infornation Center, U.S. Geological 
Survey. In order to drive the previously described energy 
balance nodel the following paraoeters would have to be 
neasured, estinated, or calculated in advance for every 
point on such a grid: 

Inconing solar radiation (Q>g) 

2. Albedo (a^ and a ) 

Q q 

34 Inconing longwave radiation 
Air tesperature (T.) 

w ■ 

5. Bind speed (u) 

64 Surface roughness length (Zq) 

7. Ihecsal conductivity (K) and vol^^tietric beat capacity 
(pC) of surface material (can vary with depth) 

8. Air pressure (F) 

S. Vapor pressure of air (e.) and at surface (ej 
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CALCOLAIZOH OF MODEL ZFFOl EABIMEIEBS OVEB A SgBFACE 


At this stage o£ sodel building# we do not sinulate 

1 

realistically all of the above facaceters. The aii here is 
to' outline the develcfient of an energy balance eodel for 
rough terrain. In particular# we sust specify# or 
interpolate froa a set of aeasureaents# albedo# wind speed# 
roughness length# and the theraal characteristics of the soil 

ovier the grid# although we calculate variations in albedo 

I 

with soil wetness and illuaination angle £163* Infocaation 
on{ albedo and perhaps roughness length are available froa 
reaote sensing# although it is necessary to sinulate incoaing 
solar radiation in order to estinate albedo. A lack of 

I ■ ■ ■ 

knlowledge of the surface thernal characteristics is not a 

I ■ 


sericus obstacle to use of the nodel at this tine# as one of 
thie innediate applications could be in inferring these 
thernal characteristics in cases where surface tenpecature is 
known from thernal inagery. Calculation of wind speeds over 
abuntainous terrain renains a sericus problen# although 
rbcent work on neso-scale wind field nodelling indicates that 
calculations which are accurate enough to be useful night be 
possitle [11, 13# 30]. 

Variation of air tenpecature with elevation is specified 

on input to the nodel, with a default specification of an 

anbient atnospheric lapse cate of ^.0065 n**^. Air 

tenpecature at altitude z can be calculated froa neasuced 

tenperature. T_ at. altitude z„; 

^ - m m 


T = T - .0065 (2-z ) (11) 

z m in 


Ihe present aodel does net consider .air drainage* but the 
coaputer prograa is so strtetured that such algocithss* if 

t 

developed* could easily be added. 

iir pressure 9 at altitude z is derived froa aeasured 
pressure at altitude z^ frea the hydrostatic eguation* 
the eguatioo of state* and eguaticn (11): 

dP/dz s -gp^ 

Pa = 

P = exp{C (ga )/(.0065B) ] ln£ (1 -.0065 (z-z J )/T ]} (12c) 

III a m an 

Ihe relative huaidity of the air and of the soil surface 
are censidered constant over the grid. Specific huaidities 
are calculated froa eapirical eguations relating saturation 
vapor pressure to teaperature £27]. 

Solar radiation is the aajor input parameter which 
varies widely over the grid* and an efficient method of 
calculating the radiation incident on every grid point is 
essential to the aodel. many slope radiation models have 
been developed £ 1h* 15* 19* 20* 24* 25# 45* 49 J and we 
present the pertinent eguations froa these only briefly. 

They are fully outlined in the original sources. 

Initially tKe beam and diffuse components of solar 
radiation on a horizontal plane at the earth's surface must 
be estimated or aeasured. Under clear shies* these can be 
calculated from the dust and water vapor coateut of the 
ataospbere and the solar geometry. The beam cadiation is: 


(12 a) 


(12b) 
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Bcooks £63 has given appcoxinations for the absocption 
and scattering coefficients: 

k, « •.089 tBfl/1013) ; »* - .174 (p^B/20) .* (14) 

fl w 

kj * -.083(fld) .« (15) 

Typical valaes for the dost and haze content d range 
froa 0 » 2 , for very clear at aospheres, to as high as 3.0 £15]. 
At present vork is underway to define vavelength'*dependent 
absorption and scattering coefficients fron aeasureaents of 
global radiation at different solar zenith angles. The 
optical air nass B is deternined by Rasten*s foraula £23}: 

B = 1/CcosZ ♦ . 15 (90-Z*3.aa5)-i -2*3 ] (16) 

where Z is the solar zenith angle in degrees. The bean 
radiaticn Q* incident upon an unshaded slope is: 

Q* = e*P£hjj+k^] (17) 

In calculating the radiaticn on an actual slope a test oust 
be made to deternine whether the sun is shaded by a local 
bofizcn. If so, Q* - 0. 2 , the solar zenith angle, and Z*, 

the solar angle aeasured fees ncraal to the slope, are given 
by; 

cos2 = sin$ sins + cos4> cos5 cosw (18) 


cosZ* s CosS cosZ «■ sinS sinZ cos (A-E) 


(19) 





Th« aziauth o£ the sun fzoa south, negativa when east of 
south# is given by: 

cosA s (sis« cosZ • sin5)/(cos« sinZ) (20) 
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Since the cosine function is even# this eguation gives tvo 
possible angles foe the aziauth of the sun# A oe •A. If ui 
is negative# then the negative value foe A is chosen; 
otheevise the positive value foe A is used. A foue-teen 
Foueiec seeies gives the sclae declination in eadians and 
the eadius vectoe to about four significant figuees: 

6, e * 2 accost 2ir (J-1) i/366] ♦ b^sin£2Tr (J-1) i/366] (21) 


Values foe the ccsine and sine coefficients foe each seeies 
aee in Table 1. 

The diffuse cadiaticn slope factoe is cos2(S/2) £24]# 
and the diffuse radiation g is the sum of the dounscatteced 
eadiaticn and backscattezed cadiation; 





g = 0.5£ 1-exp (kg) J (cos* (3/ 2) 

♦ 0.5a^£1-exp(kg) ] ♦ a^Q] 


( 22 ) 


In equation (22) ag and a^ refer to the beam and diffuse 
# albedo on the surrounding region. In cases vherc the solar 

i > . : ■ " ■ : ; 

zenith angle Z< 50®# ag= a^# but at larger aaglcs# the beam 
albedo oust be corrected for glint effects £ 16 ]. 



Calculations for horizon angles# slope# and exposure 
oust be made at each grid pcint. The method presented here 

, / : ;v D 


is diiiscent ftoa that o£ iiliiaas et al. C49]« If sore than 
one siaulation for the sane -topographic surface is nade, it 
is note efficient, because the conplete horizon and slope 
calculations are stored bn a disk or tape file, and need not 
be reconputed for different son paths. Slopes can be 
calculated front 

tans * C ( 3 Z/ 3 x)* *• (3*/3y)*]‘ * (23) 

The azisuth of the slope ficn south is given jay: 

tauE = (3 2/3 1) / (3 y) (24) 

This function is defined even if the denominator is zero, 
and if the signs of both the numerator and denominator ace 
known, the function is defined over the range to r. The 
origin for the and y-coordinates is assumed to be in the 
southwest corner of the grid, but the eguatious and the 
program can be easily modified for other orientations. 

The partial derivatives az/3X and az/ay must be 

evaluated by numerical means, and it is necessary to choose a 

method which can be adjusted when the point in guestion is 

near a boundary of the grid. He use the following group of 

nethcds: (i) if the point is on the boundary, we use a 

second-order forward difference approximatiou; (ii) if the 

point is in the row or column next to the boundary, uc use a 

seccnd'order central difference appcoximation; (iii) if the 

point is in the third row or column from the boundary, we use 

a fourth>ocdec central difference approximation; (iv) if the 

point is greater than three rows or columns from the 
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boundacy, we use a sixth-ocdet central diffeceuce 
appcoxiiatioa« Coaputational fociulas for each of these 
approxiiations can be found in £4]. 

A set of horizon angles (Z^} is generated by 
calculating the tangent of the angle froa the specified 
point to every other point in the grid whose elevation is 
greater. For: exanple the calculation froa point Ca»n] to 
Ci, j] is: 

tanS£ nn,ij] == **ij"*nin^ {h«£ (i-n) (j-n) *]]-»/* (25) 

Hithin each of a set of aziButhal increnents, for exanple 
3.75, 5» or 10 degrees r the prograa stores the naxinun value 
of the tangentw^ If b is the sector nuoiber of size n, then 



b - int{tan-irt3^hj/|i^^^^^ (26) 

The naxinun value of tanS£an 4 rl:] is selected for each b. 

Froa these the horizon angles, frcn zenithr are calculated: 

r. : ' 

. . //■ 

= ctn-MtanS£ij,b]} (27) 

This sequence results in very efficient conputer code. 
Coaparison of tangents of angles, delaying calculaticn o£ 
the actual angles until the naxiauB tangents are selected, 
ainiEizes the nunber of trigcncnetcic functions that are 
coaputed. * 

Incoiting long-wave radiation froa clear skies is 
calculated by Erutsaect*s nethcd £7]; * 

1+ = (20) 


is 


where the effective dteosphecic eBissivity 

e. » 1.2^l(e /T )»/» (29) 

' C| ,0.0 

A cecent test of Brntsaect's aethod, using data froe 
southern France, indicated that its accuracy coapare|J 
favorably with that of other sethcds £28J. For calculation 
of the net Ion g-> wave radiation , the theraal radiation view 
factor is given by Beifsnyder and Lull £41 j: 

a cos* (<2^>) (30) 


where <Z|^> is the nean horizcn angle, fros zenith. The 
diffuse radiation slope factor is cos^ (S/2) # and the net 
long**wave radiation at the site is then: 




I * I+-I+ * cos* (S/2) £ (1-V^) 1+ ♦ (7^-1) ] (31) 
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BESDLIS OF SIMOIAIZON 

1 ' ■. 

Ihe coaptttec progcai resulting froa eguations (1) 
through (31) was tested on a digitized topographic grid froa 
the Chitistone pass area (latitude 61.6<*II) in the Wrangell 
Boantains, Alaska^ and located on the McCarthy C*-! 1:63360 

quadrangle. We selected the area because a liaited set of 

\\ 

test data was available £5]. The area is about 3 ks by 3 ka 
and ranges in elevation frci about 1500 to 2100 a. The 
nearby surrounding terrain includes elevations up to 2800 a, 
and this terrain coabined vith the low sun angles produces 
iaportant horizon effects. Froa the NCXC Digital Terrain 
Tapes# we extracted a grid at 63.5 a spacing which spanned an 
area of about 11 x 10 ka. Froa these terrain data# we 
calculated the horizon vectors for the study area. The 
output presented in this paper is in the fora of ccntour 
aaps; the prograa will also draw perspective plots if 
desired. Figure 1 is a topographic aap of the study area and 
the surrounding terrain. 

For the siaulation we used the input paraaeter^^ 
listed in Table 2 below. The initial soil teaperature was 
assuaed to be equal to the air teaperature to a depth of 1 a# 
and 24 hrs of teaperature evolution were siaulated to 
establish realistic initial conditions. Values for porosity# 
thersal conductivity# and roughness length were taken froa 
field aeasureaents £5] and estimated over the rest of the 
study area froa a aap of landforas [10]- Figure 2 illustcates 
the terrain of the study acea itself and the distcibution of 
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surface aaterials. IbecBal properties for the different 
surface aaterials are given in Table 3. 

Figure 3 illustrates the distribution of incoaing solar 
radiation at 0600* 1000* 1400* and 1800 hours* The effects 
of slope* exposure* and horizon can be seen in these aaps. 

For exaaple* the point labelled A in figures 2 and 3 is on an 
ead^-facing slope and receives laxioal radiation at 0600 and 
1000. At 1400 and 1800 the intensity of the solar radiation 
at A is saall. Point B is located on a north**facing slope* 
and radiation values are lo« at all tioes* especially at 1000 
and 1400. 


> Figure 4 illustrates the distribution of net all'-wave 

■ 

radiation at 4**hr intervals iron 0200 to 2200. Unifcra 
surface esissivity is assused. The najor source of variation 
is incosing solar radiaticn* as deaonstrated by the laps for 
tines Kben the sun is not up* 0200 and 2200* 

Figure 5 shcus integrated values of incoaing solar 
radiation and net all'wave radiaticn over 24 hrs. The 
sharpest gradients occur in ccn junction with north -facing 
slopes*! as at point C (which i^^^ indicated on figure 2i. 


Figure 6 illustrates surface teapecature at 4>hr 
intervals. The effects cf different surface materials are 
cbyicus cn the 0200 and 2200 maps* where the sharp 
temperature gradients coincide with the surcaGe material 
transitions. During the day* the tempecature differences due 
to topographic effects overwhelm those due to the thermal 
properties of the soil. 


DISCUSSICH OF BESOLTS 
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lables 4 and 5 pxesent a coaparison of field sucface 
teaperature aeasaxeaents with aodel xesults for five sites, 
whose locations ace indicated on figure 2. The teapecature 
data are fcoa £5], wherein an earlier version of the cliaate 
siaulator was shown to gi/e realistic results, given the 
availability o£ site-specific input data. Two abstractions 
of the data are ccapared with the aodel. 

One is aaplitude and phase ccapacison between observed 
and siaulated values (table 4). Although the aodel 
underestinates the aaplitudes, it does rank them in the 
correct order. The phases (tine of inaxiauB teapecature) ace 
also correctly ranked. The other cosparison used is the 
difference between the 2:30 a. a. and 1:30 p. a. ’ceaperatures 
(table 5), which correspond tc the overpass tines of the Heat 
Capacity Happing Mission Satellite (HCMH) . Here the observed 

/I 

differences ake less than the diurnal aaplitudcs, because of 
slope and exposure ef iects. The acdel results coccespond 
ceascnably well to the data. 


For these liaited verification data, the aodel results 
are encouraging. However, the observed and simulated 
© teaperatures do not satch exactly, and these discrepancies 

are Frobably related to four effects. The first is an 
inadequate evaluation of sooe of the surface parameters, most 
# notably roughness length and t her sal diffusivity, both of 

which are difficult tc measure. 
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lb6 second is acdel ovecsiaplif ication. The solar 
radiation attenuation calculations in the present nodel are 
non'spectral and do not include specular cefiection iron 
adjacent terrain. Furtheriioce nean daily values for air 
tenpeiature^ relative hunidity, and wind speed are used. An 
altitude correction is applied to the air teaperatutes, but 
otherwise spatial oniforsity is assuaed, and advection is not 
considered. Therefore the nodel results are soncwhat 
over generalized. 

Thirdly# the neasured surface tenperatures nay be in 

i ■ ■ ; 

error. The nethcd that was used was to place a theraistor at 
the surface but shield it fxoa radiation. Under these 
circusstances# the direction of the error would be toward 
larger aaplitudes. ^ ' 

The fourth is the resclution cf the Digital Terrain 
Tapes. These tapes were nade fron the 1:250#000 topographic 
naps and tend to give a snocthed picture of the terrain# 
thereby reducing thersal contrast. o 

Because of these linitations# a full test of the nodel 
will have to await the availability of'theraal iaagery of 
sountainous terrain, fie suggest# however# that the results 
are sufficiently encouraging to indicated that for at least 
sone applications# the nodel is adequate. 

i The model has soae ionediate uses. One critical 
paraneter which would need to be neasured over an area before 
the surface eguilibriua temperature could be accurately 



siBttlated is the surface albedo. HoNevac« the albedo canaot 
be reliably aeasured reactely unless the incoain9 solar 
radiatios at all foints on the surface can be accurately 
siaulated. Otter nan and Fraser £31 J have successfully 
■easured albedo froa LANDSAT over horizontal terrain » but# 
foe rough terrain, naps of solar radiation, like those of 
figure 3, vould be- a necessary prereguisite to resote 
■easuresent of albedo. 

In the present application, the uind speed was 
relatively unifora over the study area. If the surface 
tesperature is to be accurately sisulated over a larger 
elevation range in sountainous terrain, the air flow and 
drainage will also have to be aodeled. Alternatively, in 
cases where thernal data are available froa aircraft or 
spacecraft, the aedel could be revised to solve for wind 
speeds 

Ihe present aodel is designed for either bare soils or 
low,; grassy vegetation. To apply it to an area of brush, 
foi’est, or cropland, where a definite canopy exists, the 
radiation attenuation by the canopy wo u have to be 

considered £4dX convective processes in the canopy 

would have to be aodelled £8]. ; 

During periods of cloudy weather, the sinulaticn of 
incoiEing short* and long*wave radiation becuaes aore 
difficult. For scae applicaticns of the model, for example 
geological prospecting ftoa thermal imagery, overcoming this 
problem is not necessary, because imagery would only be 


available on clear da;s. Hovever# if the aodel «ece used to 
aid in the siaulatioo of sone surface process'over a lon 9 ec 
tiae period/ the energy budget vouid also have to be 
calculated on cloudy days. The need for supporting 
aeasucenehts of solar radiation would increase on cloudy 

I 

days. In particular, both the beas and diffuse conponents 
would need to be seasured, sc that proper coapensaticn of 
slope and exposure effects could be Bade. If, however, 
accurate aeasureaents were available, the absorption and 
scattering coefficients could be deter ained froa the 
aeasureaents instead of froa eguations (14) and (15). 
ipplication of the aodel during partially cloudy weather 
would reguire Careful analysis of the tine distribution of 
the cloud cover. 

Incoming long**wave radiation is more difficult to 
deteraine, and for scae applicaticns of the model an 
accurate estiaaticn would be quite iaportant. For example, 
aelting snow typically has a relatively high albedo in the 
short-wave range, thereby reducing the importance of solar 
radiation and consequently increasing the importance of 
thermal radiation as an ene^ source, it is probable that 
cloud temperature estimates frcrn satellite measucemonts 

j\ 

could be used to estimate an effective atmospheric 
emissivity, and this appears to be an area for useful 
research in the future. 
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lable 2 


lapot PacaMter Valnes foe Siaulatioa 
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e 


pate 

Aix Pcessuce > 

Oust conteut of ataospbeire 
Pcecipitable water 

I 

Air teapecature 
Hr relative huaidit]f 
Dfy albedo 
net albedo 

Soil surface relative 
huaidity 

Sind velocity 

Bough ness length 

Soil 

Soil tberaal conductivity 


Ibgust 4 

835,0 Bb (at 1790 a) 

0.2 (i.e. 160 part ca*^) 

20 aa 

7.4®C 

ii 

0.81 

0.20 to 0.33 
0.15 to 0v24 
0.30 to 0.70 


2 a sec“* 

.02 to .10 a 
.08 to .48 

.87 to 4*6 Ha-4 OK-i 




, lable 3 



Tteraal Properties of Surface fiaterials (see £l 9 ure 2| 


thereal 
con duct ivitj 
H i-t ®K-I 

porosity ■ 

soil 

eater 

fraction 


albedo 

AesiduaX debri 

s 





0 to .In 

.87 

.28 

1.0 

.02 

.2 to .3 

•1 to •2a 

4.5 

• 48 

1.0 



•2 to 9 a 

.71 

.10 

1.0 



Ireas of aass 
aoueaeot 






0 to .3 a 

.92 

.30 

.7 

.05 

.24 to .33 

.3 to 9 a 

1.4 

.08 

.7 



Bedrock slopes 
vith surface 
debris 






0 to .3 a 

2.7 

25 •30 

.3 

.02 

. 15 to .20 

.3 to 9 a 

1.4 

.30 

.3 



Bora ice 

j 





0 to > 1 a 

4.6 

.40 


.10 

• 15 to . 20 
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Table 4 

Coapacisoa of Tenpecature Aaplitude aad Phase 
Aaplitttde Phase 


Sihe 

^ Observed 

Slaulated 

i 

ii 

Observed 

Simlated 

1 

17.1 

i 

14.8 

1130 

1215 

2 

16.5 

14.8 

1210 

1230 

3I ' 

1813 

16.6 

0930 

1 100 

4 

31.5 

20.1 

0930 

1115 

5 

20*5 

18. 1 

1600 

1500 



xTable S 


d 


0 


0 


0 


Site 

1 

2 

3 

4 

5 


cisoa o£ 2:30 

j 

a. and 1:30 

p.B. TeB^acatares 

2:30 a. 


* 

1:30 

p.B. 

ObsecTed 

SiBttlated 

Obsec ved 

SiBulatcd 

i 

3.f 

5.2 

16.4 

19.4 

4.7 

5.2 

18.3 

19. 5 

5.0 

5.2 

20.8 

21.4 

2.8 , 

4.9 

17.0 

18.8 

S. 0 

3.1 

19. 0 

, ^ 16.9 


oiffecence 


Site 

Observed 

Si aula ted 


1 

12.9 

14.2 


2 

13.6 

14.3 


3 

15.8 

16.2 


4 

14.2* 

13.9 


5 

14.0 

13.8 




Captions fee Flgucos 

^ 1. Contoar aap ol test acea and saccounding tercain. 

Contour interval is 100 a» 

2. Enlarged contour aap of test area at contoar interval of 
2S| a, and aap of surface aaterials. Huabered points 
identify locations iihere surface teaperature data are 
available. Points aarked by let tec are referenced in 

© ' 

text and subseguent figures. 

3. Incoaing solar radiation on August 4, froa 0600 to 1800 

^ /■•■■■■"■ 

e (soiar tine). Contour interval is 50 H a~ 2 . 


e 


4. Net all*'wave radiation on August 4, froa 0200 to 2200 
hrs. Contour interval is 50 ar^ , with supplcaentary 
contours at 25 w a~ 2 . 



Net daily solar radiation (left) and net all-wave 
radiation (right)' on Augest 4. Contour interval is 
2 HJ n-2. 


Sprface equilibciua teaperature (^C) froa 0200 to 2200 
ijis on August 42 Contour interval is 2^. 


ORIGINAL PAGE IS 
OF PdOR QUALITY 







FIGURE 2 


B*drock ilap*! wtUi wrftc* dtbf«* 


(scale for all subsequent figures) 


FIGURE 5 
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